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INTRODUCTION. 


It is the purpose of this paper to develop in outline the theory of a 
system of n ordinary linear differential equations of the first order 
containing a parameter and subject to certain boundary conditions. 
Toward this end the notation of matrices is used. For the convenience 
of the reader the paper opens with a brief review of the fundamentals 
of matrix algebra and the integration and differentiation of matrices. 
This is followed by an expository discussion of the homogeneous and 
non-homogeneous differential matrix equations of the first order. 
The major portion of the treatment is devoted, however, to the homo- 
geneous differential vector equation with a complex parameter in its 
coefficient, and to the system composed of such an equation and 
suitable boundary conditions. The solutions of the equation for large 
values of the parameter are discussed and the formal development of 
a vector of arbitrary functions into a series of solutions of the system 
is obtained. The paper closes with the proof of the convergence of 
this development under appropriate conditions, which, in the ordinary 
notation, establishes the possibility of simultaneously expanding n 
arbitrary functions in terms of the characteristic solutions of a prop-. 
erly restricted differential system of the type 


ys(a) = 2 fais(2)d + bix(x)}yx(z), 
5 fasxny x(a) + Bixyx(b)} = 0, {= 3 2,. : ne? 
k=1 


When reduced to a single equation of the nth order this includes as a 
special case the expansions obtained by Birkhoff in 1908. 





1 Much of the material preceding the proof of convergence is due to Birk- 
hoff, having been developed by him in lectures at Harvard University in the 
fall of 1920. The reorganization of this material into its present form, the 
treatment of the irregular case, and the proof of convergence are due to Langer. 

2 For other developments in this field and more complete references see the 
f = 2 ~s in the Transactions of the American Mathematical Society: 

irkho 


, On the Asymptotic Character of the Solutions of Certain Lineor 
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Section I. 
Definitions.’ 


An array of elements of the form 
! 





411 G2 Gig ....Qin 

a21 ra ae ue on 

RN os Bid gta cad ove ai Azan 
' Anil An? we ne ee ae ae Ann 





in which the number of rows equals the number of columns is called a 
square matrix, and is denoted either by (a;;) or by A. Two such 
matrices are said to be equal when, and only when, every element of 
the one is equal to the correspondingly situated element of the other. 

The sum of two matrices, (a;;) and (6;;), having the same number of 
rows and columns, is by definition the matrix (a;; + 5;;), from which 
it follows that matrix addition is both commutative and associative, 
ie, A+B=B+A,andA+(B4+0C) = (A44+8)4+C. 

The product of two n rowed matrices A and B is defined by the 
identical equation 


n 
(a;;) (bi;) = ( p» aubsi), 
k=1 
on the basis of which it is easily verified that matrix multiplication is 


both associative and distributive, i.e. 


A(BC) = (AB)C, 
A(B + C) = AB + AC. 





Differential Equations Containing a Parameter, and Boundary Value and Ex- 
pansion Problems of Ordinary Linear Differential Equations, vol. 9 (1908), p. 
219 and p. 373. 

Wilder, Expansion Problems of Ordinary Linear Differential Equations with 
Auxiliary Conditions at More than Two Points, vol. 18 (1917), p. 415. 

Hopkins, Some Convergent Developments Associated with Irregular Boundary 
Conditions, vol. 20 (1919), p. 245. 
. Hurwitz, An Expansion Theorem for a System of Linear Differential Equa- 
tions of the First Order (about to appear in vol. 22 (1921)). 

Langer, Developments Associated with a Boundary Problem not Linear in the 
Parameter (about to appear in vol. 22 (1921)). 

3 For a more ample discussion of the theory of matrices see Bécher, M., Jn- 
troduction to Higher Algebra. New York; The Macmillan Co., 1907. 
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That it is in general not commutative is a consequence of the fact that 


nN n 
> aixbx; is in general not equal to 2 b;,a;,;. The rearrangement of 
k-1 k=1 
the factors in a matrix product is, therefore, as a rule, not permissible. 


The multiplication of a scalar into a matrix has simply the effect of 
multiplying each element of the matrix by the scalar. Thus if k is a 
scalar, then kA = Ak = (ka;;). Conversely, any factor common to 
all the elements of a matrix can be factored from the matrix. 

Two special matrices must be mentioned, namely 0 = (0), the 
zero matrix, and I = (6;;), the unit matrix, where 5;; = 0 when 2 +), 
§;; = 1. These matrices satisfy respectively the relations 


AO = OA=0 
and AI = 1A =A. 


The determinant formed from the elements of a matrix without 
changing the order of the array is called the determinant of the matriz. 
The alternative notations | a;;| or | 4 | for the determinant of the 
matrix (a;;) will be used. 

Given a matrix A it follows that if | A | = 0 then there exists a 
unique solution in the z’s for each of the linear systems 


n 
2 Qik Lkjg = Sse ‘= 1, 2, -- MM, 
k=-1 


where jo ranges from 1 to n. This means that there exists a unique 
set of n? quantities 2;; such that 


nr 
2 Qik Chi = 653, 1L)= 1, me, » 0am 
k=1 


i.e. there exists a unique matrix (2;;) such that 
A: (x4;) = I. 


This matrix (2;;) is denoted by the symbol A~ and is called the inverse 
of A. From its derivation it is seen to satisfy the relation AA-= I. 

Either of the relations, AX = O or XA = O, leads, on the assump- 
tion that | A | + 0, to the conclusion X = QO, as is evident from the 
theory of the systems of linear equations to which the matrix equations 
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are equivalent. It follows from the relations AA'= I and JA = AI, 
however, that AAA = JA = AI, i.e., A(A7A — I) = O. 
Hence if A is any matrix for which | A | + 0, it is seen that 


AA —I=0,i.e., AIA = I= AA, 


In accordance with the following definitions, namely 


f A(z) ‘9 a ( f a;;(2) ir) 


dA(z) Ee | (sexte?) 
dx \ dx /’ 





and 


a matrix is seen to be integrable or differentiable if and only if this is 
true of each of its elements. It is also clear that if C is a matrix of 


dC 
constants, then — = O, while for any product 


dx 
d dB dA 
—AB=A— — B. 
dx . dx + dx 
Section II. 


The equation Y’(x) = A(x) Y(x).4 


Consider any matrix of functions Y(z) which satisfies (i.e. is a 
solution of) the equation 


(1) Y’= AY. 
In accordance with the rules for determinants we have 
Yi1 <a ie yu a] Yin ie Yin 
8S gies P eats: PI apeiniien 
dx Ee Te See ae 6 eA ee tf See eS ae ee 
Yni Ynn Yni - Ynn Yni Ynn 























4 For the general theory of matrix differential equations see Schlesinger, L., 
Vorlesungen wiber lineare Differentialgleichungen. Leipzig; B. G. Teubner, 
1908. 
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whence, substituting from the system — 


nr 

, o . 

Vi = Fa UikYki, tJ, = 1,2,...0 
, =] 


(which is equivalent to equation (1)) we obtain 




















ZAik Yk1--- Lak Y kn! EE eae Yin 
k 
d| Y| fy ) ) Yon e 2 Oak Yei...% O2k Yen 
bac’ age TEE POPP he Pa ee ee Ere | Sie Et gee a nT eee 
Yni ot ae Se WY ee Ynn Yni ae oo . oe a ee a Ynn 
DU occ cvcccccces Yin 
Y21 + bk ons Ck eo Re Yon 
. ee 2) ee 
Zz Ank Yui... 2 Ank Ykn 
k k 








d 
1.€. | Yl = aul Yl + ae] ¥| +... ann! ¥. 


Let us suppose, now (7), that the elements of A(x), Y(x), and Y’(z) 
are each continuous in an interval a S x S J, and (i), that at some 
point z of this interval | Y|=- 0. Then throughout a neighborhood 
of the point in question 

d| Y | ® 


——— = = ax, da, 


| 
| Y| k=1 


or, integrating, 


Inasmuch as the right-hand side of this equation cannot vanish if 
c + 0, it is seen that the hypothesis that | Y | + 0 for some z leads to 
the conclusion that | Y | differs from zero for all x. Thus we infer. the 
THEOREM: If Y(zx) is a matrix of functions which satisfies equation 
(1), while | Y(xo) | = 0, aS 2% 3S 5, then | Y(x)| =0,aS2xb. 
A matrix of functions Y(x) of this type which satisfies equation (1) 




















soi tt a tak 5 m oats fo . - = 
Pe ss nee ey ee Ne AS hah P 


Pg 2 loath LG ate SE ESS eS 


wi ret eA FEAL ONE eR TE Hace aN Cg Katha 


2 ohh amt Rhee 
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and whose determinant | Y(z) | #0, is called a matrix solution of 
equation (1). We assume at least one such solution to exist. 

THEOREM: If Y(z) is any (matrix) solution of equation (1), and C 
is any matrix of constants (for which | C | + 0), then Y(z), defined by 
Y(x) = Y(z)C, is also a (matrix) solution of equation (1). 


38 
Proof: We have oF =— YC 
dx dz 
dY dC 
=—C Y — 
dx + dz 
_4aY 
de 
dY 
But by hypothesis —m AY. 
da 
Hence ey =z AYC, 
dx 
dY " 
1.e. —= AJA 
1.e Tn Y 


Moreover, since the product of two matrices is derived in the same 
manner as that of two determinants, it follows that | Y| = | Y||C|. 
Hence | Y | + 0, if|C|+0 and| Y|+0. QED. 

THEOREM: Given any matrix of constants Yo for which | Yo | + 0, 
then there exists a matrix solution of equation (1) which for any 
preassigned 2x, say x = 2Xo, reduces to the matrix Yo. 

Proof: It has already been shown that when Y(z) is a matrix solu- 
tion of equation (1), then Y(x)C is also such a solution, where C may 
be any matrix of constants whose determinant is not zero. Then in 
particular C may be chosen as the matrix Y~'(a 9) Yo, whereupon it 
follows that Y(x) Y-'(zo) Yo is also a matrix solution of equation (1). 
This solution, however, obviously reduces to the matrix Yo when 
zt = Zo. Q. E. D. 

THEOREM: If Y(z) is a matrix solution of equation (1) then Y(x)C 
is the most general solution of equation (1). 
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Proof: Let Y(zx) be any solution of equation (1) whatever. Then 
the relation Y(x) = Y(ax)©® €x) defines a matrix ®(x), and substitut- 
ing from this expression in equation (1) we have 


bd YS = AY, 
dx 


d@ dY 
i.e. Y ——+—®%= AY®. 
ois a a 
dY 
But by hypothesis —= AY. 
dz 
d® 
Hence Y—+ AY@= AY, 
dz 
d® 
Le. Y — = 0. 
; dx 


d® 
Inasmuch as | Y(x) | + 0, it follows that _ O, and ®(zr) =C. 
x 
Q. E. D. 
A mere interchange of the réles played by the rows and the columns 
of the matrices involved transforms the discussion carried out thus far 
for equation (1) to the corresponding discussion for the equation 
Y’= YA. The facts in the two cases may, therefore, be established 


by precisely the same methods. 
The pair of related equations 


(1) Y’ = AY, 
(2) Zz =s ZA; 


are said to be adjoint, either being the adjoint of the other. 
THEOREM: If Y(x) is any matrix solution of equation (1), and Z(z) 
is any matrix solution of equation (2), then 


Z(x) Y(x) = C. 


Proof: From the relation 


d dZ dY 
hese die tee pea 
dx id dx ‘* dx 
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it follows, upon substituting from equations (1) and (2), that 


£ ZY = (— ZA) Y+ ZAY). 
dx 


d 
Hence rs = 0, and Z(x) Y(x) =C. Q.E.D. 
z 


ConvERSE THEOREM: If Y(x) is any matrix solution of equation 
(1), and the matrix Z(x) is defined by the relation Z(x) Y(x) = C, 
where | C | + 0, then Z(x) is a matrix solution of equation (2). 


Proof: We have @ ZY =O, 
dx 


: dZ dY 
eC. —Y+Z—=0. 
oe dz T dz ° 


In virtue of equation (1), therefore, 


dZ 
—Y+ZAY = QO, 
dx 


and it follows, since | Y | + 0, that 


dZ 
— =—ZA. 
da 


Moreover, | Z| = |¥~“|-|C|+0. Q. E. D. 


Any pair of solutions Y(2) and Z(x) of equations (1) and (2) re- 
spectively which satisfy the relation Z(x) Y(x) = I, are said to be 
associated solutions. Thus if Y(ax) is any matrix solution of equation 
(1) the associated matrix solution of equation (2) is Z(xz) = Y~(z). 

A differential equation is said to be self-adjoint when and only when 
it is identical with its adjoint after interchange of rows and columns. 
A necessary and sufficient condition that the equation (1) be self- 
adjoint is readily seen to be that a;; = — a; for all 2 and 7. 
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Section III. 
The equation Y’(x) = A(x) Y(x)+ B(x). 
THE EXIsTENCE THEOREM: Given the equation 
(3) Y’=AY+B, 
where A(x) and B(x) are matrices of continuous functions, a S z Sb, 
then there exists a unique matrix Y(x) whose elements are functions 
which are continuous together with their first derivatives, a S x S 3, 
and which satisfies equation (3) as well as the condition Y(a) = Yo, 


the matrix Yo being any prescribed matrix of constants. 
Proof: By means of the relation 


Ym(x) = Yor f { A(t) Ym_s(t) + B(t)} dt 


it is possible to define the following infinite sequence of matrices 
Y,(x), Yo(x),...¥m(x), ... which satisfy the relations 


Y;(zx) _ A(x) Yo + B(x), Y;(20) _ Yo, 
Y2(x) = A(x) Yi(x) + B(a), Y2(a%o) = Yo, 


YaiG) = 4G) Yonie) + Ble, Yale) = ¥ 
Then setting Ym(2) — Ym-i(z) = Unm(z) 
in the identical equation 
Y m(%) = Yot (¥i(x) — Yo) +(Yo2(x) — Yilx)) +... +(Ym(x) — Y m-1(2)) 
we have 


Y¥m(x) = Yot Ui(x) + Uo2(x) +... + Un(z). 


Moreover, Um (xo) = O, 


while U »' (2) 


Ym (x) — ¥ m—1 (2) 
- A(x) { Yn-1(x) sp Y m-2(x)} 
soaet A(z) U m-1(x) . 
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Denoting by the symbol << the fact that each element of the matrix 
of constants on the right represents a value as large as the largest 
numerical maximum attained by any element of the matrix on the left, 
we have, since A(x) and U;(zx) are continuous in the closed interval 
asxzZb, 


U;(2) < (k), 
A (2) <<(q). 


From this it follows that 
U;(a) = A(z) U; (x) = (Ea, uf) << (nak) 
l-1 
and integrating 
U.(2) < (kna | z= ZX | i, 


n nr 
Likewise U; = AU, = (Fa. uf) <<( 3a-kna lz — 2| ), 
l=1 l=1 


oe 2 
U3(x) << ( kn202 | a =P, 


2! 


and similarly for / = 3, 4,..., 


ee ee a aime 
U1,,(2) << (in' | x = ') <ch (= a} ) 








However the infinite series of matrices 
(3 eical) 
Zz 
0 l! 


ke (end -a —}), 





converges to the matrix 


Hence the series 


co 
l=1 


converges uniformly throughout the interval a S x S 5, and since its 
terms are matrices of continuous functions the matrix Y(x) which it 
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represents is necessarily one with elements continuous for a S zx S b. 


Now 


Yn(a) = Yor f {Yn + BO} dt 
Zo 
by definition, and since the convergence of Y»(x) to Y(x) is uniform 


lim Ym(a) = Yot ff tim { A(t) Ym-i(t) + B(t)} dt, 


m= 


namely, Y(x)= Yo+ f { A(t) Y(t) + B(#)} db. 


oa) 


Differentiating we see that the elements of Y’(x) are continuous, 
asxzx<b, while Y’(x) = A(x) Y(x) + B(x). Hence the solution 
exists as stated. 

To prove the solution unique suppose that both Y(x) and Y(zx) are 
solutions of equation (3) each reducing to the matrix Yo for x = 2p, i.e. 


Y’(a) = A(x) Y(x) + B(x), Y(ao) = Yo, 
Y'(x) = A(x) Y(x) + B(x), Y(ao) = Vo. 


Now Y(x) — Y(x) = D(z) is a matrix whose elements are continuous, 
a < xb, and which satisfies the relations 


D'(x) = A(x) D(x), D(ao) = O. 


Let any neighborhood of the point x» be chosen, say | x — 2o| S 6. 
Then, if d denotes the largest numerical maximum of any element of 
D(x) for any x of this neighborhood so that D(x) << (d), 


A(x) D(x) << (nad), 
and D’(x) << (nad). 
Hence D(x) << (nad | x — x9 | ) << (nad6). 


But for some 2x of the interval at least one element of D(x) is numeri- 
cally equal to d (by the definition of d), and for this element it follows 
that d S nadé, ie. d{1— nad} S 0. 
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However, since a and n are fixed numbers and 6 can be chosen so that 
1 

5< — this leads to a contradiction unless d= 0. It follows that 
na 


there always exists a finite interval throughout which d = 0 and hence 
throughout which D(z) = 0. Inasmuch as this implies that D(x)=0, 
a < x Sb, the solutions Y(x) and Y(zx) must be identically the same. 
Q. E. D. 

The solution of equation (3) can be easily expressed in terms of the 
solutions Y;(x) and Z,(x) of the homogeneous equations (1) and (2) 
respectively. Thus, multiplying both sides of equation (3) on the 
left by Z,(x) we have 


dY 
Zi, Ca Z,AY + ZB, 
dx 


which, in view of equation (2), can be written 


dY dZ;, 
A4,— = — — Y+Z,B, 
da dx ™ 


i.e. 


d 
+e Z,Y _ ZB. 
dx 


Integrating we obtain 
ZY=C+ f Z,(t) B(t) dt, 


while the multiplication of this equation on the left by Y,(x), the 
solution associated with Z,(x), yields 


(4) Y(2) == Y,(x) C + f Y,(a) Z(t) Bit) dt. 


This is the general solution of equation (3). In consequence any 
particular solution Y(2) may be written 


Y(x) = Y;,(zx) C +f Y,(x) Z(t) B(t) dt. 
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Subtracting this from equation (4), however, we obtain 
Y(x) = Y(x) + Y,(x) {C — C}, 


from which it is seen that if Y(x) is any particular solution of equation 
(3) the general solution is given by 


(5) : Y(x) = Y(x) + Y,(zx) C. 


If in particular C is taken as C = O in formula (4), the solution which 
is characterized by the fact that Y(a) = O is obtained. This is called 
the principal solution at x =a. Since the choice of a as a limit of 
integration is unrestricted, the expression for the principal solution at 
any chosen point is at hand. 


SecTIon IV. 
The differential system 


¥"(x)+ = A(x) Y(x)- + B(x)s 
W,Y(a)+ + W,Y(b)+ = O. 


A matrix in which any row (or column) is precisely like every other 
row (or column) is called a vector. That a particular matrix is a vector 
is indicated by means of a dot suitably placed in relation to the letter 
designating the matrix in question, the dot preceding in case it is a 
vector in which the rows are the same and succeeding in case it is a 
vector of identical columns. 





A= (a;) nn oe ee ee ee 
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That -AB = -C, and that AB- = C-, are facts readily established by 
direct reference to the rules for multiplication. A matrix (a) all of 
whose elements are the same is written -A-. From the preceding 
statement it is clear that -AB- = -C-. 

Given W, and W;, any two constant matrices for which | W, | + 0 
and |W,|=+ 0, then the two corresponding homogeneous vector 
differential systems 


6) Y’(x)- = A(x) Y(z)- 
W.Y(a)- +W.,Y(b)- = O 
and 
-Z(a) W.+ -Z(b) W= O 


are said to be adjoint. 


THEOREM: The number of linearly independent solutions of system 
(6) is always equal to the number of linearly independent solutions of 
system (7). 

Proof: It has been shown that if Y(x) is any: matrix solution of the 
differential equation (1), the most general solution is Y(x) = Y(zx)C. 
From this, C = Y—(x) Y(x), and if the solution Y(x) is a vector 
Y(a)- then C will be a vector C-. 

The general solution of the differential equation in system (6) is, 
therefore, 


and since the substitution of this in the boundary conditions gives 
(8) W, Y(a)C-+W,Y(b)C: = 90, 


Y(a)- is seen to be a solution of system (6) when and only when C:- 
is a solution of the equation (8). Moreover, it is readily seen that a 
necessary and sufficient condition that a set of solutions of system (6) 
be linearly independent is that the corresponding solutions of equation 
-(8) be independent. Setting W, Y(a) + Wz» Y(b) = (p;;), equation 
(8) can be written in the form 


n 
Pik Ck = Q, = a ae 


(9) 


k=l 
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The number of linearly independent solutions of the linear algebraic 
system (9), (and hence of equation (8)) is, however, precisely equal to 
the difference between the number of equations, n, and the rank of the 
determinant | ;;|. Thus when | W, Y(a) + W, Y(b) | = 0, and is 
of rank (n — k), then there are precisely & linearly independent solu- 
tions C;-, C2-,...Cx-, of equation (8), and correspondingly just k 
linearly independent solutions Y,(x)- = Y(x)C;- i= 1, 2,...k, 
of system (6). Conversely if there are k independent solutions of 
system (6) the rank of the determinant is (n — k). 

Suppose then that system (6) is known to have just k linearly inde- 
pendent solutions. This means that the determinant vanishes and is 
of rank (n—k), and since | X,U X;2 | is of the same rank as | U | when 
| X,| = 0, and | X.| + 0,° it follows that 


| Z(a) Wo { Wa ¥(a) + Ws Y(b) } Z(b) W| = 0 


and is of rank (n — k), Z(x) being any matrix solution of equation 
(2). If in particular Z(x) is chosen as the solution associated with 
Y(zx) this reduces to the statement that | Z(a) W.+ Z(b) W.“| = 0, 
and is of rank (n — k), which implies that system (7) also has just 
k linearly independent solutions. Thus the theorem is proved. 

A system of the type (6) or (7) is said to be either compatible or in- 
compatible according as it does or does not admit of a solution not 
identically O. Its compatibility is said to be k-fold when the number 
of its linearly independent solutions is k. 

Consider the non-homogeneous system . 


Y"(x)- = A(x) Y(x)- + Bla): 


(10) W, Y(a)-: +W, Y(6)- = 0. 


THEOREM: A necessary and sufficient condition that system (10) 
has a unique solution is that the corresponding homogeneous system 
(6) is incompatible. 

Proof: It was shown (formula (5)) that the general solution of equa- 
tion (3) is 


Y(x) = ¥,(x) C + Y(2) 





5 Cf. Bécher, loc. cit., pp. 77-79. 
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where Y(zx) is any particular solution and Y,(zx) is a matrix solution 
of the homogeneous equation (1). Hence the general solution of the 
differential equation in (10) is 


Y(x)- = Y,(x)C- + Y 2)- 


The substitution of Y(a)- in the boundary conditions shows it to be 
also a solution of the system (10) provided only that C:- satisfies 


W.Y(a) +W,Y(b) +{W.Y,(a) +W, Y,(b)} C- = O. 


But this relation can be solved for C-, and uniquely determines C- 
when and only when |W, Y;(a) + W, Y,(6)| +0, that is when 
system (6) is incompatible. Q. E. D. 

Assuming then that system (6) is incompatible it is possible to 
obtain by the following procedure a solution of the equations (10) 
which is symmetrical with respect to the ends of the given interval 
aszczb. 

Let Yi(x)- and Y2(x)- be any pair of solutions of the differential 
equation 


(11) Y'(x)- = A(x) Y(x)- +3 B(a)-. 
Then clearly their sum Y,(x)- + Y2(x)- = Y(zx)- is a solution of equa- 


tion (10). But by formula (4), applied to equation (11), a particular 
choice of Y;(z)- and Y2(x)- is seen to be 


Vile) = 4 f Yale) 4 BO-a 
Ya(2)- = 4 f Yate) 2 BO-at 
b 


from which it follows that 


Fz) = 44 f Yale) Za(t) BO-at +f Ya(z) Z(t) BQ)-at 
a b 


is a particular solution of the differential equation (10). 
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Defining G(a, t) by the relations 
EB 5 Y,(x) Z,(t) whent <2 
G(z, t) = 


—4 Y,(x) Z,(t) whent > 2, 


we may write 


b 
Y(x)- = |G. t) B(t)- dt. 


In accordance with formula (5), therefore, the general solution of 
the differential equation (10) is given by the equation 


b 
Vin)<: i G(x, t) B(t): dt + Y,(2)C- . 


The substitution of this form in the boundary conditions yields the 
equation 


b 
[WG + W.G0,0} BO-at + (7. Yala) + W¥,0)}C- = 0 


for the constant vector C-. Multiplying by the inverse of the matrix 
A = {W. Y,(a) + W» Y,(6)}, 


(A-! exists since system (6) is incompatible) we see that 
: 
C-=— a+ f {WG (a, t) + WiG (b, t)} B(t)-dt. 


It follows that the general solution of system (10) is, in terms of a 
matrix G(x, t) which is defined by the formula 
(12) G(a, t) a G(a, t) we Y,(x) A" {Wa G(a, t) + Wy G(b, t)}, 


given by the expression 


b 
(13) Y(z) = fo t) B(t)- dt. 
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The matrix of functions G(z, t) is known as the Green’s function for 
the homogeneous system (6). By a precisely similar method the 
Green’s function H(z, t) for the related system (7) may be derived, the 
general solution of the non-homogeneous system 


Za) = — -Z(x) A(x) + - D(z) 


(14) -Z(a) Wo + -Z(b) Wx= O 


being given in terms of H(z, t) by the formula 
b 
(15) -Z(x) = f- D(x) H(z, t) dt. 
a 


THEOREM: If G (2, t) and H(z, t) are the Green’s functions for sys- 
tems (6) and (7) respectively, then 


G(z,t) + H(t,7) =0, t4z. 


Proof:§ Let B- and -D be any two vectors of the types indicated 
- and consider the two systems 


Y’-= AY-+B.,, W., Y(a)-+ W, Y(b)-= 0, 
‘Z = —-ZA+-D, ‘Z(a) We+ -Z(b) Wy'= O. 


Multiplying the differential equations respectively by -Z on the left 
and Y- on the right and adding we obtain 


‘27 + PY: @ BB + DY-. 


an equality which upon integration yields 
b b 
(16) ‘LY: x = a -Z(t) B(t)-dt + f -D(x) Y(x)-dz. 


Now from the boundary conditions we have 
-Z(b) = — -Z(a)W.' Ws Y(b)- = —W. WwW, Y(a)-, 
whence 
-Z(b) Y(b)- = a) Y(a)-, 
i.e. ‘“ar-| = 0. 





6 The proof by direct computation is not difficult though somewhat laborious. 
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b 
But -Z(t) = f -D(x) H(t, x) dx 
and Y(x)-= f G (zx, t) B(t)-dt 


from formulas (13) and (15). 
Hence equation (16) is equivalent to 


b 6b 


(17) ; : f -D(x) {H (t, x) + G(a, t)} B(t)-dx dt = O. 


This result, having been obtained without reference to the nature of 
B- and -D must, moreover, hold for all possible choices of these 
vectors. We shall proceed to choose a particular set apposite to the 
proof in hand. 

Let any point (20, to) of the region a S x,t S b,x +t, be arbitrarily 
chosen and draw the surrounding small rectangle As whose sides are 
t=t)+ At and x=2)=+ Az. Then choose B- and -D so that 
b, (t) =0,whenl +7, 6;,(t) -OinAs,  5;,(t) = 0 outside of As, 
d;(x) = 0, whenk + io, d;,(x)+O0inAs,  d;,(x) = 0 outside of As. 


For this choice equation (17) is equivalent to 


w+At wmtdr 
(18) ¥ di,(x) { higig(t, x) T Jioio(@s t) } b;,(t) dxdt = 0, 
t-At z0-Ar 


and inasmuch as d;,(x) b;,(t) + 0 in As, clearly {hijo(t, 2) + Gioio(2, t)} 
changes sign, i.e. vanishes somewhere in this region. Now let Az and 
At approach zero. Then in the limit we have 

hinig(to, ao) + Jing Xo, to) = 0. 
From this it is seen that G(xo, to) + H(to, x0) = O, and since (2p, to) 
was any point not on the diagonal it follows that 


G(x, t) + H(t,z2) =O, xt. Q. E. D. 


By direct reference to formula (12) it is readily seen that the Green’s 
function possesses the following characteristics: 
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i) The elements of G(z, t) are continuous in 2x except for z = t. 
Along this line (x = t) there is a jump of unit magnitude in the 
elements of the principal diagonal, i.e. 


G(z,x — 0) — Giz, 2+ 0) = I. 


ii) For any given t, G(z, t) satisfies equation (6) in x, except along 
the line x = t, i.e. 


OG(a, t 
nn = A(x) G(z,t), rt. 
Ox 
iii) For any given t, G(a, t) satisfies the boundary conditions of 
system (6) in 2, i.e. 


W,G(a, t) + WG, t) = O. 


Conversely we have the 


THEOREM: The dependence of G(x, ¢) upon the variables x, t is com- 
pletely determined by the characteristics (i) to (iii) above. 

Proof: Suppose @ (2, t) isa matrix possessing the characteristics (i), 
(ii), and (iii). By (i), (ii), J(2, t) defined by J(z, t) =G (2, t) — G(z, 2) 
is continuous for all z,a < x < b. Moreover, J(z, t) satisfies (iii) and 
is therefore a solution of system (6) inz. But system (6) is incompati- 
ble by hypothesis. Hence J(z, t) = O 


and G(a, t) = G(z, t). Q. E. D. 


It is readily verified that a further set of three characteristics which 
completely determines the dependence G(x, ¢) upon the variables may 
be obtained by interchanging x and ¢ and replacing system (6) by 
system (7) in the discussion above. 

The fact established by this theorem should be carefully noted. 
While the choice of the pair of associated solutions Y; and Z;, on the 
right-hand side of equation (12) is not unique, yet the entire function 
G(a, t) is independent of that choice. 
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SECTION V. 
The formal solutions of the equation 
Y’'(x)+ = { A(x)A+ B(x)} Y(x)-. 


Returning to the homogeneous equation let us consider the nature of 
the solutions when the matrix coefficient of Y- is made to depend 
linearly upon a parameter \ which is free to take on all values in the 
finite complex \ plane. To this end we shall study the equation 


(19) Y'(x)- = { A(x) d + Ba) } Y()-, 


where A(x) and B(x) are matrices of continuous functions, by making 
the assumption that the equation has a formal solution 


-. Mawar 


l 1 
Po(x)- + x P(x): -+- \2 P2(x)- +... 


where @ is any chosen constant. 

If n = 1, equation (19) can be directly integrated and is seen to have 
an actual solution of the form (20). The passage to formulas (45) can, 
therefore, be made directly in this case, and hence we shall assume 
in the intervening work that n 2 2. 


z 
Setting f y(x) dx = T(x), and substituting the form (20) in the 


equation (19), we obtain the formal identity 
eae 1 ) 
y(x)e™t (2) Po(x) . -- r P(x) ‘ + ear + ” aie (2) P)(x) + eee 


] 
= {A(x)A + B(x)} } Pae)-+ - P,(z)-+... , 
from which it is seen, upon equating the coefficients of \, that 
(x) Po(a)- =A (x) Po(a)-. 


This is satisfied by a vector Po(x)- not all of whose elements vanish, 
when and only when 


(21) | asj(x) — 5;; y(a) | = 0, 
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i.e. when 7(x) satisfies equation (21). For any given 2, say z = 2, 
however, the left-hand side of (21) is a polynomial of degree n in y(2»). 
Hence the equation is satisfied by n roots y1(%0), Y2(ao),...¥n(2o). 
We shall assume that for the case in hand these roots can be grouped 
into the nm functions (x), ye(x),..-.yn(x), which satisfy the three 
conditions 


(i) (2) continuous, 72=1,2,...n, aszxxsb, 
(22) (ii) y;(x)  y:(x) for 7 + 2, 
(ii) yz) +0, «= 1, 2,...n. 


Clearly the last condition can be fulfilled only if | A | = 0; we shall 
assume this to be the case. 

Consider now a change of the dependent variable in equation (19). 
Setting Y(x)- = &(x) Y(x)-, where &(x) is a matrix whose elements are 
continuous as well as their first derivatives, a S x S b, and |@| +0, 
the equation becomes ®’Y-+ 6Y-’= {An + B} OY-, that is 
Y’-={@" Ab + BS — S"G'} Y- It is evident here that if the 
coefficient of \ is not originally J then no such change of variable has 
the effect of reducing it to 7. In the subsequent work it is desirable 
to transform A(z) into the matrix R(x) given by R(x) = (6;; y;(x)). 
By such a transformation the formal solutions (20) are carried into 
others with the same functions 7;(x) in the exponents. 

Since the ®(z) in question must satisfy -! Ab = R,’ it must fulfill 
the conditions 


(i) Ab = OR, 
(ii) | (x) | + 0. 


To satisfy condition (i) the elements of ® must be solutions of the 
algebraic system 

nr 

Dik Oki = Oi Vn 4 J, = 1, 2,...0, 
while the possibility of fulfilling this condition by means of a matrix 
(x), none of whose columns contain only vanishing elements, depends 
upon the existence of a solution for each of the n linear systems 





7 For fuller discussion see, for instance, Bécher, loc. cit. Chap. X XI. 
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n 
(23) > las — dix vi,.} oxi, = 9, 
k=l 


obtained by giving to jo successively the values 1, 2,...2. This 
means that each of the determinants | aj, — 5:x Yj, |, jo = 1, 2,...m 
must vanish, and inasmuch as the functions 7;(2) were chosen as the 
roots of equation (22), which is precisely the condition that the de- 
terminants in question do vanish, the existence of a matrix ®(z) of the 
type desired is established. 

Condition (ii) requires further, however, that this matrix ® be one 
for which | &(z) | + 0. 

Suppose | @(a)| = 0. Then the system of linear equations 
nr 


> viz 0% = 0 admits of solution by means of a set of quantities 1, 
k-1 


Vo,. . . Un, at least one of which differs from zero,” i.e. there exists a vector 
V-# O, such that @V-= 0. But then ASV -= O, and in view of the 
relation (i) it follows further that PBRV-=O. By continued repeti- 
tion of this reasoning it can be shown likewise that ®R*V- = O for all 
values of k. Consequently ®{ col + ¢,R+c2R?+...+ en-1 me) V-=0. 
But on the hypothesis that y; = yi when 7 + 2 it is well known that 
the determinant 








ly vv Y1 

ly. ye ...¥2"" 

A ey ee te ae = 0. 
ae ein 


Hence the system of equations 
cote + ey? +...+enay” = ki 


Cot CiY 12 + COV n? +... Cas y= Ras 


udmits of a solution in the c’s not all zero for any choice of the quanti- 
ties k; not all zero. Some choice of the set ¢o, ¢1. . .¢n-1, can therefore 
always be made to satisfy the relation 


feol a C} R +. _ -} Cn—1 | a = (6;; k;). 


8 Cf. Bécher, loc. cit., pp. 47. 
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Then &(6;;k;) V- = 0, 


nr 
le. 2 gy k, vx, = 0, for all choices of the set ki, ko,...kn. But by 
l=1 


construction some 2, say 2;, is not zero. If then the k’s are chosen so 
that k; = 0, 7 + jo, k;, £ 0, it follows that ¢;; = O forall 7. Inasmuch 
as no column of ® contains only vanishing elements this involves a 
contradiction. Hence the hypothesis | © (x) | = 0 is not tenable, and 
the ® in question fulfills condition (ii), i.e. | &(x) | + 0. 

By a change of variable, then, equation (19) can be given the form 
Y’-= {RX + B} Y-, where the matrix B, being given by the relation 
B =" B — 4’, is a matrix of continuous functions. Supposing 
this to have been done we may drop the dashes and consider the equa- 
tion in the form 


(24) Y'(x)- = {R(w)vA + B(x)} Y(a)-. 


Since we have found that the Y(x)- of (20) may be any one of the n 
vectors Y;,(2)- obtained by replacing I'(2) by I';(2) the further dis- 
cussion might well be carried through for each of these vectors indi- 
vidually. However, if the matrix E(x) is defined by the relation 


E(x) = (6;; 1), 


it is readily seen that the matrices P;(x) can be chosen so that the j“ 
column of the matrix 


(25) Y(x) = | Pot + : P,(z) +... E(x) 


is precisely the general column of the vector Y;(x)- Hence all cases 
are simultaneously treated by the consideration of those formal solu- 
tions of the matrix equation 


(26) Y’(x)= (R(x)\ + B(z)} ¥(z) 


which have the form (25). 
Inasmuch as E’(x) = (dy;(x) 6;; eT) = \ R(x) E(x), the formal 
substitution of (25) in (26) gives the identity 
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1 1 ee 
(27) rot rth Pet... b RE + | Rt Pt. b Em 


jar + BI } Pot Pi... | E. 


Equating the coefficients of \ we obtain the relation 


PoR = RP», 
i.e. pp Vi = Vi py? , 
whence it is seen that py = 0, whenj +7. Equating the coefficients 
of \” we obtain in similar manner the equation 
P:R + Po’ = RP+ BP», 
1.€. 
(1 (0 (1 0 

(28) Die V3 + Di = Vs Dig + be bir Diy 
whence it is seen, upon setting 2 = 7 that 

nr 

0) ‘(0 0 

DS} os = bik Pip = bi; Dij' 


J b;;dz 
namely, that p‘j = f; e* 


Having obtained in this manner all the elements of the matrix Po 
we find further from the equation (28) for 2 + 7 that 
(0) 
ay _ uP, 


(29) sane 
Yi~—vYi 


; l 
Moreover, equating the coefficients of 1 we see that 


P.R + P, = RP.+ BP,, 


‘ n 
2 1)? 1 
ps} {yi- v3} + ps; iat * bix py ’ 


whence it follows upon setting 7 = 7 that 


(30) Ps; - = bik Diy : 
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Since the other quantities under the sign of summation in this formula 
have all been determined, p‘;’ may be found by means of a quadra- 
ture. When this has been done the matrix P; has been completely 
determined. 

It should be observed that we have thus far required only that the 
matrices R(x) and B(x) be continuous. ‘To determine the elements of 
P, however, we have first, for 7 = 7, the equation 


nr 

1), 1) 

as py + 2 bi Phy 
Py = k=1 





b 


Yi Yi 


and since this formula implies the existence of p{}’, whereas it is seen 


from (29) that ps} has a derivative only if this is true of the matrices 
R(x) and B(x), we cannot proceed to the determination of P2 if R(x) 
and B(x) are merely continuous. 

Let us suppose then that R(x) and B(x) both possess continuous 
derivatives up to and including those of order k 2 1, but that perhaps 
one of these matrices possesses no such derivative of order (k + 1). 
If, in particular R(x) and B(z) possess infinitely many derivatives we 
may take k = o. The derivatives of P; up to and including that of 
order k are now seen to exist from formulas (29) and (30). 


l 
Equating formally the coefficients of iin the identity (27) we 


pl 
have 
PUR + Pa = RP, + BP,-1: 
1.e. 
(u) ‘ baie) 
~1)" ae 
Dif lyj— yj=- Pi er bik Dei» 
saciiee (u-1) 
—p; +2 by pig’ 
whence pry = Soh for 7 + j, 





<2 Peete 


; ; > I 
and again, equating the coefficients of — we have 
a 


Py R+ P, = RP, + BP, 
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1.€. pyr? ly;— vi} + py” == bsx pe, 


whence it is seen, upon setting 2 = 7 that 


nr 
pi’ = = dp}, 

k=1 
namely that pif ) can be determined by means of a quadrature if the 
quantities on the right are known. 

The determination of the elements of p;;’ by means of these formulas 
depends, therefore, only upon a knowledge of the elements of P,_; 
and upon the existence of P.. . Moreover, it is seen that in general 
P,, possesses one less derivative than P,-:. Inasmuch as P, has 
already been determined, and was seen to possess k derivatives it is 
clear that the matrices P, for u = 0,1,..., (¢ +1), may be suc- 
cessively determined, and that in general the matrix P;,; is merely 
continuous. 

If & is finite, we can, therefore, determine a differentiable matrix 
S(x), given by 


(u) 


| Po(x) + . , Pal) +...+ 5 Pale) re 


as (ane bp 4 at > +t) 


which will satisfy the equation 


(31) S(x) 





S’(x) = {AR + B} S(z) +3 ie — B(x) P;(x)} E(x), 
i.e. an equation of the form 
(32) S’'(a) = jane + B(x) + -, L(x, d) S(x), 
where L(x, \) is a matrix each element of which is rational in \, 


: 1 
with coefficients continuous in 2, given by power series in (: : 


If on the other hand k = ~, as — terms of the infinite series 


(2) 


a . 


pi = pg + + 4+. 


















LE EE tay Cpa 
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as are desired may be used. These series are not in general con- 
vergent. Nevertheless the formal matrix 


S(x) 


l 
Po(x) +y Pyi(a) +..... ta) 


(1) 
(ors Ee + Ya | 


which is found in this case, formally satisfies the equation (26). 


With the convention that the term in xi is to be omitted if k = 


formula (32) holds in all cases. S(x) will be called a formal matrix 
solution of equation (26) regardless of whether k is, for the case in hand, 
finite or infinite. 

It should be observed that each pet ) is not wholly determined but 
contains a single arbitrary constant of integration, independent of z. 
This arbitrariness corresponds to the fact that any convergent power 


series ¢; in (*) with constant coefficients may be multiplied into each 


column of S(x) without thereby destroying its property of being a 
formal solution in the sense (32). 
Using the notation [Wo]x, or [yo], for an expression of the form 


Vi y y(2,r) 
Yot+ “Gites cs Ga 





where y is bounded for | \| large, we have 


Fits 
S(x) = (| ce" . |, )Ee. 


It is clear that an alternative form is 


(33) S(x) ‘aad («, oh" 7'2)+ Ba) CAP ), 


where B(x) = f b;; dz. 


By similar considerations the equation Z’ = — Z{A\ + B} may 
be transformed and a formal matrix solution 7'(2) for the resulting 
equation 
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(34) Z'(x) = — Z(x) (R(x)r’9 + B(z)} 
be obtained. This 7(x) has the form 
(35) T(z) = (6 e-APi(2)-Bite) (5, .) 
and satisfies an equation of the form 
T’(x) = — T(x) } ARC) + B(z) +3 Mle, d) . 

Moreover each row is a formal solution of the vector equation 
(36) -Z' (x) = — -Z(x) {R(x) d+ B)}. 

Considering differentiation as merely a formal process defined by 


the usual rules it is readily seen that the differentiation of the formal 
matrices S(z) and 7(x) is permissible. Hence we have 


d 





a, US = TS’ + T'S 
= r\aR+ B+ Lh S— r}R+B+4ut S 
1 

3 TiL—M}S. 

Since 
TS = (3 Ci tala lSsale ¢; et a Lda, ) 
sy oi) og) 

= (é o; el Pi@)-Tale) }+ Bi@)- Bi) | $3; + — +....+ ee t), 

we find upon differentiating, and removing the exponential factor, 





C; ¢; INyi—s) + bj — bis} cn Ur 2 | 


+e + {fm 
» eee hah 
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1 nm 
fad lOsale ns — Me} [Sali 


Equating the coefficients of XX we have __ 
6:0; { (y¥i— Vs) FF + (bij— Bis) 5:;} = 0, 
from which it follows that of; = 0 when i+ j. Again, equating the 


1 
coefficients of . we have 
(y;— Yi) of + (bj;— bii) of + of; Or = = 0, 


from which it is seen, upon setting 2 = 7 that o})' = = (), namely that 
(1) 


o}; = constant. The relation shows on the other hand that of = 0 
when 1+). Equating to zero successively the coefficient of each 
1 
individual power of > it is found in the same way that o{) = 0 when 
$+ 3, of) = = constant for uw = 1, 2,...,4¢ —1). It follows that 
o? (k-1) 
A{Ty-P;}4+B;-B; o; o4;(2, d) 
en (rnin ea 
= (6; ¢; 54; (1x1) + 5 E" (0;;) E, 

where the coefficients of [1],_1 are constants. 

Now for any choice of the set of series c; it is clearly possible to 
choose a set Z; such that ¢;¢;[1],1 = 1, 7 = 1, 2,...n. The formal 


solutions S(z) and 7(x) corresponding respectively to these values of 
c; and ¢; are closely related. They are called associated formal solu- 


1 
tions and satisfy the relation T(x) S(x) = I+ ye EO (o;;) E. 


Since in particular c; may be chosen as c; = 1 it is seen that there 
exists a formal solution having the form 


(37) S(x) = (eTI@+Bj@ [6,,)). 


In accordance with the definition above the associated solution is 
given by the set of ¢;'s which satisfy the relations ¢;= 1. It is 
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apparent, therefore, that the solution associated with (37) is of the 
form 


(38) T(x) = (e>Fi@-Bi® [5,)). 


Section VI. 
The relation of the formal solutions to the actual solutions. 


It was observed in the preceding section that the formal matrix 
S(a) either satisfies equation (26) only in an approximate sense (i.e. 
satisfies (32)) or, if k = © satisfies it only formally, since the elements 
of S(x) are in that case infinite series which are not in general con- 
vergent. S(x) is, therefore, not a matrix solution of equation (26), and 
its significance requires further investigation. 

Consider the actual matrix which is derived from S(x) by retaining 
in the latter only the first (m + 1) terms of its elements, where m is any 
positive integer not exceeding k. This matrix S(x) may be written, 
in accordance with formula (37), 


S(x) = (eATj(2)+Bj(2) [5sj)m) 
. = PG, d) EG), 


(ei [5:j1m), 


where P(x) 
and is seen to be analytic ind. Since | e7/6;; | + 0, it follows that 


|S|+0for|\’|> N. If the formal solution S(x) is written as 
P(x, X) E(x), it is apparent that we have the formulas 


S(z) = 1 ~ ai (Q] ph sc, 





S(x) = 1 + ri (Q] P| S(x), 


from which we obtain upon differentiating and substituting from the 
equation (32), the relation 
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ar 1 1 1 fe 

S’(x) = }r- oa pt har+B+41| — ym ((2l Py 
11+ a 1 P45) 


It is apparent from this that S(z) is a solution of a homogeneous 
differential matrix equation of the type 


(39) $'(2) = bar +B+ = ®(x, d) S(z). 


l 
The elements of ®(z, \) are power series in (*) with coefficients 


which are continuous in 2, and are seen to be convergent, since these 
elements are rational in X. 
But equation (26) can be written in the form 


(40) Y"(a)= } Ria) + B(z) +00, »f ¥(x) — a P(e d) ¥(z), 


and considering this as a non-homogeneous equation we know, in 
virtue of the developments of page 62, that its solutions, i.e. the solu- 
tions of (26), are given by 


. | ieee 1 
- C cise : 
Y(x) = S(xz) C + f S(x) T(t) am &(t, d) ri dt, 
l.e. by 
“): temteien 4 f Giz) FU) HE, ») Y@ dt 
_ 


where 7(x) = S-(x) and where the lower limit of integration, which 
has been omitted, may be chosen at pleasure for each of as many parts 
of the integrand as desired. 

Substituting for Y(a) in equation (41) its equivalent as given by 
the form 


(42) Y(x) = (x) Sr), 


we have further 
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(43) U(x) = S(x) CT(2) — = f S(x) F(t) &(t, ») U(t)S(t) F(a) dt. 


THEOREM: If the functions y;(x), 7 = 1, 2,.. .n, satisfy the relations 
arg {y;(x) — y:(x)} = hij, 2,7, = 1,2,...n, where each h;; is a con- 
stant, then there corresponds to each sector bounded by two adjacent 
rays R{d{y;(x) — yi(x)}} = 0° a choice of the lower limits of integra- 
tion which is such that for \ within the sector and | \| > N, and for 
any continuous matrix U(t), each element of the matrix 


¥(@, d) = f S(x) T(t) S(t, ) Ut) S(t) F(x) at 


is less numerically than KM, where M is the largest numerical maxi- 
mum attained by any element of U. 
Proof: Writing 


S(x) = (s:(x)) E(x), T(x) = E-(z) (t:;(z)), 


we have 


He) = ( .. tate Un abies 


h,l, p, @, r=l 


Sor(t) eM Tr) —Tr(z)} t,;(x) it) 


n z 


as rS { v,(&) — v,(€) }de 
3 aa 8 i wt? (t, a, d) mS 


where, for || > N,|o | < «@"M for all i and j, k“” being a 
positive constant. 
Consider a particular sector and any element 


Zz. 
> rf {4 —-7,@ }de 
he é : yt) dt. 


If R{A{-y,(€) — 7,(€)} S 0 for \ within the sector and any ¢ it is so 
for all €, and, provided that ¢t S z, the integral 





8 The notation R {¢} is used to indicate “ the real part of ¢.” 
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Zz 
 r>S {y, © —7,© } ae 
4 é . wo) dt 


a 
is less numerically thank“? M. However, if R{A{y,(£) — y,(£)}}> 0 
for \ within the sector and any &, it is so for all £, and, provided that 
t > x, the integral 


Zz Zz 
rS Ly, © —7,(€) Jae 
f “its wo dt 
b 
is similarly bounded. Consequenily if x,,(A) is defined by the rela- 
tions 


Khp = @ if Ridtyna) = vr(é)} } S 0, *,, = 6 otherwise, 


the numerical value of each element of the matrix 


Zz Zz 
2" (Sly @-7,@ ae 
¥(x, r) = ( ie J. . wo"? dt 
h Kh 


,r=l 


is clearly less than KM, where K = (b — a) = k“”, for X within the 


h,r=1 
sector and |\| > N. Q. E. D. 
Assuming the limits chosen in the manner above, a ¥(z, A) corre- 
sponding to each sector and to each U is uniquely determined. More- 
over, we have 


(44) U(x) = B(x) CT (2) — = V(x, ). 


Consider now the particular solution, Yo(x), of equation (26) which 
satisfies the relation Yo(a) = S(a). Since S(z) as well as the coefficients 
of equation (26) are analytic in \, Yo(zx) is likewise analytic in X. 
Moreover we know from page 56 that every solution of the equation is 
of the form Y(x) = Yo(2) D, where the elements of D are constants 
with respect to x. Substituting this form of Y(z) into equation (41) 
and fixing z and solving for C we obtain the relation 





C= T2)} Yo(xo) + xn [Se T(t) P(t, \) Yo(t) dt¢ D, 
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n 
Cij = 2 pin dn;, 
hl 


where the quantities p,;, are analytic in X. 

Inasmuch as Y(2) is a non-identically zero solution if D + O there 
will exist such a solution for which C = O provided the determinant 
|p, | vanishes. If on the other hand this determinant does not 
vanish then there exists a solution Y(x) = O corresponding to every 
C + 0. : 

Suppose the determinant | p;,| = 0. Then there corresponds to the 
choice C = O a solution Y(x) for which the matrix U(x) defined by 
Y(x) = U(x) S(x) satisfies the relation 


pe 1/*% : Ni vp, (&)—-vr(€) } ae (hr) 
oe * (E. Je ie it) 
Kher 


But we know that for some x, say 2, and for some 2, 7, say 20, Jo, 
Us (to) = M. Then since 


x Zo 
° Mf {-,(@)-v,(€) } de 
|Z foe mT" a? (wo, t) dt| << KM, 


h, r=1 


Kh 
KM 
we have for the 20, 79 element, M < \m? 
K 
1.e. u(1-X)<o. 


Inasmuch as M is positive and \ can be taken arbitrarily large this 
involves a contradiction, and proves the hypothesis untenable. 
Hence | p;,{=- 0. Accordingly there corresponds to every choice of 
C + O, and hence to the particular choice C = J, a unique D + O, and 
hence a solution Y(x). By equation (44) the U(x) for this Y(x) must 
satisfy 


U(x) = I ae W(x, r). 


We infer from this the inequality, 


KM 
M<1 Tm? 
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whence it follows that for X sufficiently large, M is less that 2 and 
¥(x, X) is uniformly bounded, for ) in the sector in question and further 
|X| >N. Accordingly we have 








Y(z) = ( 6:3 + tae ») (s:;(z)) E(2), 
i.e. Y(zx) = ( 8i;(a) +— = ») E(z). 


Since the matrices Yo(x) and D are analytic in \ the same is true also 
of Y(x). It is seen, moreover, that | Y|+0. Hence, when the 
functions 7;(x) satisfy the conditions of the theorem on page 83 there 
exists in every sector of the type described above a matrix solution 
Y(x) whose elements are continuous in 2 and analytic in \, and are 





the same as those of S(x) to terms in But by construction the 


,=" ° 


mS 1 
elements of S(x) are the same as those of S(x) to terms in \m* Hence 


we have the 

THEOREM: Given any formal solution S(z) of equation (26) in which 
R(x) and B(x) both possess derivatives up to and including those of 
order k, then in each sector of the complex plane within which none of 
the quantities R{d{-y,(x) — y,(x)}} change sign there exists an actual 
matrix solution Y(z) which is continuous with its first derivative in z 
and analytic ind, and is throughout the sector identical with S(zx) to 


terms in =. 
N 


In virtue of this there exists in any sector of the type described a 
pair of actual associated matrix solutions of the form 


¥ (x) = (e¥j(#)+8;(2) [6,,],_1) 
(45) Z (x) = (e->Fs(@)-Bs(2) [5,],1) , 


provided only that the matrices R(x) and B(x) are differentiable 
k terms. In all cases the functions [6;;],_,; which occur in the expres- 
sions for Y(a) and Z(a) are, of course, respectively identical, to terms 
in 1/** with the series [6;;] in the formulas (37) and (38). 
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In the course of the deduction of forms (45) it was assumed through- 
out that n 2 2. Direct integration of the equation, however, shows 
that in the case n = 1 there exist solutions which are of this form over 
the entire plane. Accordingly formulas (45) may be used in every 
case. The explicit forms 


Y (x) - (onsoreaye 5; .s on d) ), 


Z(x) = (-arste-2ie hi + aa d) | 


with 9;;, ¥;; bounded, hold if R(x), B(x) are differentiable. 








Section VII. 


The characteristic values of the system 


Y’(x)-={ R(x)A+ B(x)} Y(x)- 
W.Y(a): + WsY(b)-= 0. 


It was shown on page 65 that a necessary and sufficient condition 
that the vector system 


y’. = {RA+B} Y: 


(46) W, Y(a)- +W, Y(b)- =0 (|W.|+0, |W,| +09), 


has a solution is that 
(47) |W. Y(a) + W, Y(b)| = 0, 


Y(x) being any matrix solution of equation (26). Let us make the 
specific form of this condition in any sector of the \ plane apparent by 
substituting a Y(x) which is analytic in A for | \| > N, and which has 
within the sector the form 


(48) Y¥(x) = (e1i*%; [5,)), 


as determined in the preceding sections. Choosing the a of formula 
(20) as a = a it follows that Y(a) = ([(6,;)). 

Substituting from (48) in the determinant on the left of equation 
(47), (call it D(A)) we have 
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D(A) = | Wa ([5id]) + We (55 [5;,) |, 
1.e. 


(49) Dr) = | (w] + [7] APeOrBO| | 


In this as in subsequent formulas when the expressions in question 
are determinants the letters r and c rather than 72 and j are used to 
indicate row and column respectively. 

Due to the fact that those and only those values of \ which satisfy 
equation (47) are characteristic values, i.e. values which yield solu- 
tions of system (46), equation (47) is known as the characteristic 
equation of the system in question. __ . 

The introduction at this point of quantities 5; and 85; will do much 
toward simplifying the further discussion. These quantities are de- 
fined by the relations 

= 6;; when R{AT,(b)} S 0 =0 when R{AT,(b)} SO 
543 oj 
= (0 when R{AT;(b)} > 0, ! = 6;; when R{dAT;(b)} > 0. 


It is to be noted that 5;; and 6;; are functions of arg \ alone. 

Suppose now that the rays® R{AT,(b)} = 0, i = 1, 2,...n, are 
drawn in the plane of the parameter \ and let one of the sectors 
bounded by a pair of. adjacent rays of this kind, RIAL 2(6)} = 0, and 


R{AT, (b)} = 0, be denoted by oy. Then if 6 and 6%? are respec- 

tively the 555 and 6;* for some (any) particular value of \ within ¢y, 
* x 

we have the identities 6;; = 6%), 6;; = 6°%?, for all \ in the interior of 

the sector in question. It becomes apparent in virtue of the relations 


[w’®?] when R{AT,(b)} <0 
(a) (b)] are sali 
(50) [whel + [wie] erPe) +B) = [ao] edPe() + Bel) 


when R{AT,(b)} > 0, 


and the fact that R{AT;(b)} + 0, i = 1, 2,...n, that D(A) takes, for 
any A in the interior of such a sector the form 





9 A half-line issuing from the origin of the \ plane will be referred to as a ray. 
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“ * b * 
(51) D(A) = | [whe] bee + [wie] bce erF eH) + Bet) |, 


Factoring from the determinant the exponential factors, any one of 
which occurs in each element of an entire column, the further alterna- 
tive form 


n OK 
(52) Di) = I ol APE (O) + Be O) }5xk | [ww] 52, + [ro] ot | 
kel 
is obtained. It becomes necessary at this point to differentiate be- 
tween certain types of conditions which inherently characterize any 
particular system of type (46). 

The conditions of the system will be said to be regular (i) if n = 1, 
or if, when n 2 2, arg I';(b) + arg {+ I';(b)} when j + 7, and (ii) the 
boundary conditions are such that each of the determinants 


_ ** 
kl kl b kl 
(53) WH = | wi BE + wO 5h? 


differs from zero. 

The conditions of the system will be called zrregular if either of these 
conditions is not satisfied. In the further discussion we shall consider 
first the case of regular conditions and then the particular type of 
irregular case which results from dropping the regularity condition (i) 
above. 

Case I. Regular Conditions. In this case the rays R{AT;(b)} = 0, 
1= 1,2, ...n, are distinct and divide the plane into 2n sectors of 
type %,;. Let us fix the attention upon any one of these rays, say 
‘the ray R{AT,(b)} = 0, and let 7,, and o,, denote the abutting sectors, 
the former being that within which R{AT,(b)} <0. For X on the 
ray R{AT,(b)} = 0 the v™ column of D(X) consists of terms whose ~ 
form cannot be abbreviated by means of relations (50). Retaining 
therefore the original expressions for the elements of this column we 
have 

* OK 


Dr) = | fw] 8% + [w®] {5% + 5,,} dBc + Bel) | , 


For \ either within ¢,, or within ¢,,, on the other hand, the form of 
D(A) is given by formula (51), the difference in the value of D(A) in 
these sectors being accounted for by the difference in the values of 
5.- and ne i. 2 ee 
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But it is readily seen that 
(rr) _ (ur) (Pr) _ (us) 
(54) Sec = See — Ser, dice = See + dev. 


Hence if S is any sector which includes the ray R{AT,(b)} = 0 and 
overlaps a part of each of the sectors 7,, and ¢,, we have, for a d in S, 
the following — 


D(r) = | [w'?] 5” an [w°”] g(u) eFe(b) + Belb)| for \ within ¢,,, 


Dr) = | [w@I) 5 (wn + [w°?] he ; + 5,,}erFe) + Bcl)| for X on the 
ray R{dT,(b)} = 0, 


D(r) aid | [w®] 5m) Ser} oh [w?] {5 + 6, ,} erTc(b) + Be(b)| for X 
within ¢,,, 
It is readily seen from this that D(A) is given for any \ in S by the 


formula 


(55) D(r) poe | [w'?} a$ (ur) abe fw {5 4. Sey }erPelo) + Bc()| 


2K 
Factoring from this determinant the product I el ATK) + BECO) } BG, 


k=l 

we have 

(56) D(d) = Terlra +BL0} 94) Dory, 
k=1 

where 


DO) = | [wo] a + (wo L3E + 8,4} edPo+ BHO) | 


Now D(a) is seen to have elements consisting of a single term in every 
column but the »”, the elements of that column being binomial. Con- 
sequently the expression of the determinant as the sum of two others 
is possible, i.e. 
are ” 
Dd = | [v5] 2 + [wo] 9” | + | ol tae? — — ber} + 
fw If on + Sep} |edF oO) + Br) 


or, in view of relations (53) and (54), 


(57) Dd) = (W"] + [W™] PrO+B-0, 
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Since the roots of D(X) = 0 and those of D(A) = 0 are the same we 
have, therefore, as the characteristic equation 


(58) (we) se (wen) orl u(b)+ Bu() _ 0. 


This yields, since W“” + 0, the equation 


we ns we 


ATy(b)+By(b) __ ee Pe 
. ¢ | a lid 





where ¢ is here introduced as a generic symbol for functions which 
approach the limit zero uniformly as | \ | increases beyond limit. 

Solving for \, and observing that log (K + e) = log K + «, we see, 
therefore, that every characteristic value which lies in sector S is of 
the form 





— yw) 
(59) A, — B,(b) + log yon) +e€(A,) + 2pm ; 


a 1 
~ T(b) 
where p is a positive integer. Moreover, since Y(x) is analytic in A 
throughout the sector S, || > N, the same is readily seen to be true 
of «€ also. 

Consider now a small circle of fixed radius r drawn about the point 








ae we 
3 B,(b) +~ log (vr) +2 prt , 


r,(b) WW 


for any given p. Then for a proper choice of origin in the \ plane 
(see page 98) this circle lies entirely within the sector ¢,, or ¢,,, and 


a <rfor||>N. Also if p is sufficiently large the point 








1 —w” 
r.) B,(b) + log yen +e€e-+2 pm 


lies within the circle for all \ on the circumference. Consequently 





bn ine B,(b) +1 algerie {| 
arg oe T, (b) a 0g wer ° inte 


increases by 27 as \ describes this circumference, and just one root of 
equation (59) is accordingly seen to lie within the circle. Since it is 
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readily verified that every such root is a characteristic value it is seen 
that (59) determines such a value for every p which is sufficiently 
large, i.e. that for large values of \ the characteristic values lie ap- 
proximately along a line parallel to the ray R{AT,(b6)} = 0, the dis- 
tance between two adjacent ones approaching as a limit the finite 
length 27/|T,(b) |, as |X| increases indefinitely. From the deriva- 
tion of this result it is seen moreover, that a similar sequence of 
characteristic values lies near each ray R{dAT,(b)} = 0, and that no 
further distribution of characteristic values exists. 

Case II. A Type of Irregular Conditions.!° Let us suppose now 
that the regularity condition (ii) is fulfilled, that n 2 2, but that 


(60) arg T',(b) = arg {+T'(b)} 


for the pair of valuesi = »;, 7 = ve. We have then a case of irregular 
conditions, and while we shall restrict the discussion to the case when 
the relation (60) holds for only a single set of values 7, 7, the reasoning 
to be employed is typical and may be applied with equal success to the 
cases in which a greater number of the points I’;(b) are collinear with 
the origin of the complex plane. It is only for the sake of brevity that 
the simplest, rather than the most general case which results from 
dropping the regularity condition (i) is treated. 

A review of the discussion applied to the case of regular conditions 
readily shows that the methods employed there apply equally well to 
. the case in hand and yield the same results in any sector of the \ plane 
which does not contain the line R{AT, (6)} = 0. It is, therefore, 


necessary to consider here only the distribution of characteristic 
values in a sector S containing this line. Along the line in question 
neither the expressions for the elements of the pth nor of the yi column 
of D(X) can be contracted by means of the relations (50). Accord- 
ingly it is found that D(A) takes the form 


* KK 
(61) D(A) = | [we] 822 + [wre] 1562 + Ser, + Ser,} eT + Bel) | 
* Aok 
throughout the entire sector S, 6% and 6% representing the quantities 
5: and dee for.A on the ray R{dT,,(6)} = O which liesin S. Factoring 





10 For the discussion of a differential system representing a different type of 
irregular conditions see Hopkins, loc. cit. 
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from (61) the product It el 14 (0)+ Bi) } 3) we have 
k=l 


n -$ar.(b) +B,() 50) — 
(62) D(A) = TI e' — ke Dy), 


where 
ohOK 


D(A) =| [w®Js¢ © + [w®] {69 + 6,, etn © +3n®) 
+ 5,,, ertoa®) +B ral) by 
_ CV, ° 


The determinant D(A) is seen, therefore, to have monomial elements 
in each column except the v{* and the vi", the elements in these 
columns being binomials. It is clear, therefore, that D(X) can be 
expressed as the sum of four determinants each containing only 
monomial elements, i.e. 


(63) D(A) = | [w?] 6 ” + [w°?] § ®) el + Raine — 30} 


+ [w 115 + Sey | orl (0)+ By, (0) os 
| Moho ~ be,,} + [wo It 30 + Ben,} [ers @ Bra) 4 
I wa — bag Sev,} 


ro Ha? + Sev, + be», ler! Py (0) +1 v9(0) +B, (01+B,,(0) ° 


Let us again denote by ¢a,, and o,, the sectors abutting on the ray 
R{AT, (6)} = 0 in which RiP, (6)} < 0 and > 0 respectively. It is 
necessary to consider the case in which arg r, (6) = arg I’,,(b) and 
that in which arg I’, (b) = arg {—T,,()}. 

Sub-case A. arg I’, (b) = arg T,,(0). 

In this case the quantities Rivr,, (b)} and R{AL,, (6)} are of the same 
sign throughout the sector S, and the relations 


oa 


as 
(oe + Ber, + der, = dee = de” 


2K OK 


(vr) (v) (uv) 
| bce Sev, es Sev, = bce = See 


(64) 
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are readily verified. Defining the determinants W; and Wz by the 
formulas 


W, (a) Rt ee Sey (b) ay Sev, 
65) = | wre | d+ wre + de} | 


W .= | w? ser + Sev, } + w? (ge) _ o, } | 


we have D(X) = D(A) where, in view of (64), D(A) satisfies the equa- 
tion (56) and is given by 


(66) D(r) ca [wo] as (W,] ed Pn(0)+By(0) 4 WV) ohlv2()+B,2(0) 
-}. (wer) e{ Bo1(b) +12) }+B,,(b)+ By (0) 


Sub-case B. arg T,,(b) = arg {— Yr, (b)}. 
In this case RiaQr,,(b)} has throughout S the sign opposite to that of 
RiP, (b) }, and it is found that 


: 5S” +8, = bf, = 88” + 8... 
(6 ) +*k KK ek 
bce” — Ser, = See = 5o6” — Sey, - 
Accordingly we have 
D(r) ie (W.] ie [W]e Fv1O)+ By, (0) re [W/o ph Pva(0)+ By 20) 
+ [W]e Py (6)+T (6) }+B,,(b)+B,,() ; 
If now we define D(A) by the relation 
D(r) oP vq(0)-By2(b) sa D(d) 


it is readily verified on the basis of formulas (67) that D(A) again 
satisfies equation (56) while it is given in this case by 


68) DA) = (W] + [Wy OOF 4 [BJM nO-Ba 
of (Wer ert Pv: (b)-T,4(b) }+B,,(6)—B,,(0) , 


Let us suppose now that the notation has been so chosen that 
r,,(6) 
lr, (5) 
equations (66) and (68) through by [W™”], (recall that W™ + 0 by 





|T,,(6) | 2 |T.,(6) |, and set 








=P, Then upon dividing 
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hypothesis) it is seen that the characteristic equation, i.e. D(A) = 0, 
is in each case of the form 


(69) 1 + feJet A. [ee] et Ty, (0) a" [es] el itr par», (0) iat! 0, 


where ¢;, 2 = 1, 2, 3, are complex constants, cz + 0 by hypothesis, and 
ris a real constant r = 1. 

Wilder ** has shown that the roots of this equation are asymptotically 
represented by those of the equation 


(70) 1 4 C} orl (0) +. Co eh FO) +4. Cs elitr JAT,, (0) = 0, 


and has discussed the distribution of the roots of this equation. We 
shall proceed to this discussion, observing, however, first that when 
| T,,(6) | and | (6) | are commensurable a far simpler treatment is 


possible. In that case r = . where p and gq are integers and the 
AT,, () 


(a 


equation (70) is an algebraic equation of degree (p+ q)ine @ 
Accordingly it has (p + q) roots, i.e. 


e @ =a j= 1,2,...,(p+q) 
from which it follows that 
a: <a : : 
(71) Az = r,,(b) {log a; + 2kwi} . 


In this case, therefore, the characteristic values which lie in sector S 
for | | > N are asymptotically represented by a set of points which are 





2qr 
spaced at intervals of length iT (| (p + q) lines (not neces- 


sarily distinct) parallel to the line RT, (b)} = 0. 

When I,,(6) and I,,(b) are incommensurable no such simple treat- 
ment is possible. The distribution of the characteristic values may 
be obtained, however, by Wilder’s procedure,!* which follows. 

Setting AT,,(b) = 2 = x + ty we have as the equation (70) 


f(z) = 1 + ere? + ese* + egel*"1* = 0, 


11 Wilder, loc. cit., p. 423. 
12 Cf. Wilder, loc. cit., pp. 420-422. 
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and it is readily verified that there corresponds to each choice of an 
arbitrarily small positive constant x, some value of x, say x = X, 
such that 


(i) |1 —f(—z)| <x, 
(72) : f(2) 


(ii) le3 — teade! <x, fora SX. 
F é 








We shall assume that x is chosen sufficiently small to preclude the 
vanishing of f(z) outside or on the boundary of the region | z| S X. 

Now f(z) is analytic throughout the entire finite plane, and hence 
it is possible to find in any interval of the Y axis, however small, some 
point y = yo which is such that the line y = yo contains no zero of f(z). 
Let y = Y; and y = Y2 be any two such lines and consider the rect- 
angle K bounded by them and the lines x = X andxz = — X. We 
shall determine the number of zeros of f(z) within K by observing the 
increase in arg f(z) as z describes the perimeter. 


We have arg f(z) = sin“! af fl2)j where I{ f(z)} denotes the coeffi- 


| f(z) | 


cient of V—1 in the expression for f(z). Moreover, for y = constant 


I{ f(z)} has the form 
It f(z)} = d,e* + dee + de elitr} . 


where the coefficients, d;,2 = 1, 2, 3 are real constants. The finite 


Tif(z)} 
| f(z) | 





zeros of I{ f(z)}, and hence those of being the roots of the 


equation 


d; +- doel?-1}2 + dse"” = (), 


are, however, separated by the finite zeros of the derivative of the left- 
hand member, namely by the roots of the equation 


do{r — 1} + dzre* = 0. 


Since this equation is satisfied by at most one value of z it follows that 
Ii f(z)} 
| f(z) | 





vanishes at most twice, and consequently that arg f(z) changes 
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by less than 37 as x varies between the limits x = —X andz = +X 
along a line y = constant. 

Because of the relation (72i), however, we know that for every 2 
on line x = —X, f(z) lies within a circle of radius x about the point 
z = 1, and hence that as z moves along this line arg f(z) changes by 
less than 2 sin-! x. Similarly relation (72ii) shows that arg{f(z)/e"*”*} 
changes by less than 2 sin-! x as z moves along the linex = X. From 
the identity 

arg f(z) = arg ell*"}2 + arg aoe 
it follows, therefore, that the increase in arg f(z) as z moves along the 
line x = X from y = Y; toy = Y¢ lies between (1 +1) {Y2— Yi} + 
2sin— x and (1 +r) {Y2— Yi} —2sin“ x. Consequently the in- 
crease in arg f(z) as z describes the perimeter of K lies between 


(itr) {¥Y.—Yi} + 67+4 sin x and (1++r) { Y2— Yi} —67—4 sin-x, 


and accordingly the number of zeros located in the interior of K must 
lie between 
(1+r) {¥2—Yi} 
2a 





ae 
+3+ —- sin x 
T 
and 


(+r) {Ye—Y¥i} _ 
2m 


Since x can be chosen arbitrarily small, however, this means that the 
number of characteristic values between any two lines y = c, andy = 





ae 
3—— sin x. 
T 





: 1 r 
¢:+ lis at least 1 — 3, and cannot on the other hand exceed 


l+pr 
2a 


Summarizing the results it is seen, therefore, that the characteristic 
values in sector S lie in a strip bounded by two parallels to the line 
R{AL, (b)} = 0, and that they are so distributed throughout this strip 


that for || > N no more than three lie between any two lines which 





t+ 3. 





2 
are at a distance d < ; ~— from each other and are perpendicular 


to the line R{AT,:(b)} = 0. 
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Section VIII. 
The formal expansion of an arbitrary vector. 


It was shown in the preceding section that both under regular 
conditions and under the type of irregular conditions discussed the 
characteristic values for system (46) are numerable and cluster about 
the point\A = ©. Denoting these values by Ay, Ao, As,. . .it is possible, 
therefore, to assign the subscripts in such manner that | Am | S | Amy: ]. 


Moreover then lim | |Am | = ©. 
m= 


Assuming that system (46) is simply compatible at the characteristic 
values there exists for each of these values just one solution of the 
system in question and just one solution of its adjoint system. These 
solutions for \ = A, will be designated by Y“(x)- and -Z“(x) respec- 
tively. 

Now if \ = 0 is a characteristic value for system (46) let the para- 
meter be changed by setting \ = \ + ¢, c being a constant. Equation 
(46) then becomes 


Y’(x)- = { Riz) + B(z)} Y(a)-, 
where B(x) = c R(x) + B(2). 


The characteristic values of the system thus modified are \ = \,—C 
and it is clearly always possible to choose ¢ so that \ = 0 is not a 
characteristic value. No loss of generality is entailed, therefore, by 
the assumption, which will be made, that \; + 0 for any k. 

Writing the equation (46) for \ = \,; in the form 


y"(a)+ = B(x) YO). + R(x) Y™ (2): 


and. considering this as a non-homogeneous equation we have from 
page 67 


< b 
(73) Y (2). =r, f G(x, t) R(t) Y(t)-dt, 


where G(z, t) is the Green’s function for the system 
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Y’(x)- = B(x) Y(zx)- 


(74) W. Y(a): +W, Y(b)- = 0.8 


In precisely similar manner we have from the adjoint system 
b 
Z(@) = —d f Z(t) R(t) H(a, t) dt, 
a 


or, substituting from the relation 


G(t,z) = — H(a, t). 


b 
(75) Z(t) =, f -Z (x) R(x) G(2, t) dex. 


a 


Consider, now, the integral 
b »b 
J= f f -Z (x) R(x) G(a, t) R(t) Y (t)- dx dt. 
In view of relation (73) we have 
b 
An = f -Z (x) R(x) Y (x) «dz, 
while it is seen from (75) that 
b 
, J = i Z(t) R(t) Y (8) - dt. 


By subtraction, then 


(An— dz) J = 0, 
and it follows that J = O provided k+ 1. But then A.J = 0, ie. 
b 
(76) Y -Z (2) R(x) Y (x) -da = O, for k + 1. 








13 Since \ = 0 is not a characteristic value of system (46) system (74) is 


incompatible and the Green’s function exists. 
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It can, moreover, be easily shown (see page 107) that when the 
system (46) is simply compatible at every characteristic value, this is 
not true for k=l, i.e. 

b 
f-2@ R(a) Y?(x)-da + O, for any l. 
a 

Let us suppose now that an arbitrarily chosen vector F(x)- can be 

developed into a series of the form 


(77) F(z): = ¥ cy Y(a)- . 
kel 


Multiplying both sides of this equation on the left by the vector 
-Z (x) R(x) and integrating term by term we have formally 


b b 
f -2@ R(x) F(a)-dx = Fes| -Z (x2) R(x) Y (x) -dz, 
a k=l a 


which in view of relation (76) reduces to 
b 


b 
(78) f -Z (x) R(x) F(x)-dx = ¢; f -Z (x) R(x) Y(2x)-dex. 


Inasmuch as the matrix on each side of this equation is one all of whose 
elements are identical the equation may equally well be written 


» 
ii zh (x) ya(x) f(x) dx (1) = 
, 
: J 2h) (2) Yn (2) yh (2) de (1), 


whence 


b 
f = xh (x) a(x) fa(x) dx 
hel 


; , 
[ 22@) n@ W@ dz 
v’ h-l 





Ci = 
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Consequently we have 
b ) 


; 
J Zn (x) va(x) fala) dx 
(79) F(z)-= r oe Ge) 
J 22 @ n@) WO az 

















If, therefore, F(2)- may be developed into a series of form (74) which 
converges in such a manner as to legitimatize the processes above, 
then (79) is a necessary form for the series in question. 

Thus far it has been stipulated only that Y“(x)- and -Z“ (zx) be 
respectively solutions of system (46) and its adjoint for \ = Ax. But 
each of these systems is homogeneous, and if Y“"’(x)- and -Z“”(zx) are 
any particular solutions then c Y(x)- and¢-Z”(z) are also solutions. 
Having chosen a definite pair Y“(x)- and -Z“*(x) we have then 


. 


[pn 

J 228 @ n@ he a 
F(x)-= r 7 , P(g). 

Jz = ezh” (x) vala) yh” (x) de 


\ @ 














If in particular ¢ is chosen so that 


if bah vn(x) yh’ (x) dx = 1, 


and the vector ¢-Z”, for this value of ¢ is associated with Y“- so that 
the choice of one implies the choice of the other, we have, on dropping 
the bars over the letters, 


) 

co b 

(0) F(z)-= L fz E al” (x) yq(a) fa(x) dx ) Y(a)-. 
k=1 


In using this formula it must be remembered that -Z™ (x) is deter- 
mined as soon as the particular Y‘”’(x)- is chosen. 








102 BIRKHOFF AND LANGER. 


In a precisely analogous manner it may be found that a development 
of the type 


F(z) = LG -Z (2), 
k=1 


which converges in such manner that it may be integrated term by 
term after being multiplied by any of the vectors R(x) Y(x)-, must 
necessarily coincide with the expansion 


b 
co 


FQ =E | J LAG n@ We ade} -2C) 


Inasmuch as the matrix G(z, t) is not a vector the methods outlined 
above do not apply directly to the problem of expanding the Green’s 
function. We shall proceed, therefore, as follows. 

Let G,(a, t)- denote the vector each of whose columns is the j“ col- 
umnof G(z, t). From the relation (80), F(x) being replaced by G;(z,t)-, 
it follows that for any value of t, G;(z, t)- can be formally developed, 
the series obtained being 


b 
(81) G,(x, t)- = Z. | > zh (2) vn(x) hj (x, t) dx Y(z)-, 
kml | Ani 


In view of equation (75) written in the form 


(25 (t)) = ry | > zh (x) va(x) gn; (2, t) de. 


=| 


a 


(81) reduces to 


co _(k) 
ee 
G,(2, t)- i ” z(t) Y)(x)-, 
ae 


from which we obtain the relation 


wo , (k) (k) 
(82) (2, =O () ys (2) 


k=l Nx 
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But since -Z‘”(#) and Y“(z)- are vectors 
3 B 
Y™(z)--Z() = ( > yi?(a) so) = n (yf (x) 2} (2)). 
l= 


Hence we have from (81) 





co ylk)y.\.. 7k) 
(83) ie bo FS 


k=l Nk 


Having obtained this formal development of the Green’s function 
it is possible to state the 

THEOREM: If the development (83) converges to G(z, t) in such a 
manner that a uniformly convergent series is obtained by multiplying 
(83) on the right by any matrix of continuous functions and integrating 
term by term, then any vector F(x)- the elements of which are con- 
tinuous and have continuous first derivatives, and which satisfies the 
boundary conditions 


W F(a): + WsF(b)- = 


is represented by a convergent development of the type (80). 
Proof: The relation 


F’(x)- = B(x) F(a): + Ce): 


defines the vector C(x)-, and inasmuch as F(x)- is then a solution of 
the non-homogeneous system 


Y'(x)- = B(x) Y(x)- + C(z)> 
W.Y(a)- + W,Y(b)- = 


it is given by the formula 
F(x): = { G(x, t) C(t)-dt (see page 20). 


Substituting for G(z, t) its series development we have 





io y ()(y). -Z (t) CH: 
rey = SE eer . 
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or, integrating term by term, 
b 
oo 1 
(84) F(z)} = FY“). § — f Z(t) C(t) «dt 
k=l NA k ne 


Now given any matrix of type D-, and another of type -A-, we have 


n 


D--A+ = (3 d; a) = (nad;), 1.€. 
l=1 


(85) D--A- =naD- 


Hence (84) reduces to 


\ 


a 
oo f 2 2(*)(t) e,(t) dt 














F(z). = 2 3 » ¥(2)-, 
k=l 
\ Nk 4 
the series on the right converging uniformly to F(2)- Q. E. D. 


In formula (12) the explicit expression for the Green’s function 
G(x, t) for an incompatible system of type (6) is given, Y; and Z; being 
respectively a solution of the differential matrix equation Y’(x) = 
A(x) Y(x), and of its adjoint. Now for any X, not a characteristic 
value, system (46) being incompatible is precisely of type (6). Conse- 
quently we have from (12) 


(86) G(a,t,d) = G(a,t,r) — Y(x, NA7AA){W.G(a, t, +) + 
W 2G(b, t, d}, 
for any A not a characteristic value. 

Let us investigate now the functional dependence of G(z, t, 4) upon 
the parameter \, choosing as the solutions Y(a, A) and Z(t, A) a pair 
which are analytic in X. Recalling that when the determinant of a 
matrix A is denoted by A then A-!= “3 , where A;; is the cofactor 


of the ij element of A, the explicit form of A-1(A) is seen to be 


A-"(X) sia (7x2), 


14 This notation holds when n = 1 provided it is agreed that when A is a 
matrix of one element then A,, = 1 
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Now it was seen (page 87) that at a characteristic value D(A) 
vanishes. Let us assume further that the characteristic values of 
system (46) are all simple so that D(A) vanishes at each of these values 
to only the first order. Due to the relation 


(87) | Di | = D™, 


which is familiar from the theory of determinants,’® it follows, then, 
that D;;(A,) cannot vanish for all 2 and 7, since in the alternative case 
the left-hand side of (87) would vanish to the order n, while the right- 
hand side vanishes only to the order (n — 1). Hence (Dj:(Ax)) + O, 
and A-(A) has a simple pole at each of the characteristic values. 
Moreover, system (46) is simply compatible at the characteristic 
values as was assumed above. It is easily verified (i), that the poles 
of A-(A) actually persist in G(a, t, \), and (ii), that G(z, t, \) has no 
other singularities. 
Denoting by G™ (z, t) the residue of G(z, t, \) at X = Ay we obtain 
from (86) by familiar methods 
G0 = YO (a) (Dix) {(WaYO(a) — WYP O20 


soos 








AwmAp 


d 
which, inasmuch as — Do + 0, is of the form 
dy may 


(88) G™ (x, t) ma y™ (2) om Z(t). 
From this it follows that for fixed t, G“ (a, t) is a solution of the system 


(89) bs. = {R(x) \,+ B(x)} ¥(x) 
W.Y(a) + WiY(b) = O. 


But if Y(z) is a particular solution of the differential equation of this 
system every solution is of the form 
Y(x) = Y(zx) C. 


Substituting in the boundary conditions as on page 64 and denoting 
by (d;;) the A which results from Y = Y we have 








15 Cf. Bécher, loc. cit., p. 33. 
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(90) (d;;) (ci) = O, 

or 2 dj cy; = 0, 1,), = i Sa. 
l=1 


Now it was shown that (D;;(Ax)) + O, i.e. that the rank of D(A,) is 
(n — 1). It follows that the solution of the system of linear equations 


(91) x dc; = 0 ¢=1,2,...” 
l=1 


is unique, namely that the ¢,; of system (90) is the c; of system (91) 
for all 7. Consequently the solution of system (89) is unique and of 
the type Y(x)C-. One solution of (89) is known, however, namely 
Y(z)-. Inasmuch as G“(z, #) was also seen to be a solution in z it 
follows that each column of G(x, t) must be the same as the general 
column of Y (x). except possibly for a factor independent of zx. 


Accordingly gs; (x, {) = c(t) y\” (x), 


1 
from which G” (a, t) = -— Y(z)--0 (2). 
n 


But from (88) it is also seen that G(z, #) is, for fixed x, a solution 
in t of the system adjoint to (89). Since this solution is again unique 
and -Z‘”(#) is a solution we have 


g(a, t) = y(a) c(t) = cb (a) 2 (0), 


i.e. H(t) = o® 2 (4), 


Hence gi; (x, t) = & y{(x) 2°) 

or 
c*) 

(92) Ge, t) = — Y~)-- 2, 
nN 


the value of c“) being as yet undetermined. 


Writing equation (46) in the form 
Ya). = {R(x) + Bix)} Y(x)-+ {rn— A} R(x) Y(2)-, 


and considering this as a non-homogeneous equation, we have 
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b 
Y(x)-= — {ry — dg} f G(x, t, s) R(t) Y(t) - dt. 


Inasmuch as 
lim {X — Ax} G(z, t,) = GM (z, 2), 


AmAp 


this yields , 
Y(z). = — f G(x, t) R(t) Y(t) -dt, 


and substituting from (92) we see that 


b 
© i(k) 
Y(z). = — | — y® (2)..Z 4) R(t) Y(t) dt, 
7 


a 


But in view of the relation (85), 


b 
y(z)- f -Z™(#) R(t) Y™(t)-dt = 


h=1 


zz 
n | > 2h (t) valt) yh (dt? Y“(2). 
a 
It follows that 


1= —¢ J 2\”(t) y,(t) y(t) dt, 
=) 
a 


or upon associating with Y“(x)- the proper -Z“ (x) (see page 101) 
that c = -1. 
Hence 


(93) G” (2,1) = — i yx). -Z (2). 
7 


This result enables us to deduce an expression for the sum of m 
terms of the series development of an arbitrary vector F(x)-. Thus 
we have from (77) and (78), (-Z’ and Y‘~ being associated in the 
manner stated) 
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b 
ri -Z (a) R(x) F(x)-dx = ¢;(1). 

Multiplying this siiciines by Y‘(x)- on the left it becomes 

b 

oY (a)-(1) = f Y(a)--Z(t) RW) F(t)-dt, 
which in view of (85) and an reduces to 
b 

(94) ¢(Y(2). = — } G (x, t) R(t) F(t)-dt. 

a 


This expression for the J” term of the formal development was de- 
duced upon the hypothesis that ; is a simple characteristic value. 
If at \; a number of characteristic values Xj, d;,,...Az,. coincide, we 
shall define the term of the formal series which corresponds to this 
value of A to be 


b 
= J G(x, #) R(t) F(t) -dt. 


In every case then we have a formal series which is completely de- 
termined and of which the / term is given by formula (94). It is to 
be noted that G (a, t) for the case in which }, is not a simple charac- 
teristic value is not given by the left-hand side of (94). It must be 
computed for every such case individually. 

Now let p; be any contour in the \ plane which surrounds ); but no 
other characteristic values. Then by use of the relation 


1 
ae f G(z, t,) dd = G(z, t) 
2712 > 


we have from (94) 
b 
] o 
¢ YO (a): = — — f f G(x, t, ) dd R(t) F(t) -dt, 
“aT * 
a Pp} 


or, upon interchanging the order of integration and choosing a contour 
C .» which encloses the characteristic values \;, \2,. . .Am and no others, 


a 
m 1 
(95) yc Y(z). = — + Fh oe. t, ») R(t) F(t)-dt dy. 
bot 271 Sie 
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SECTION IX. 
The Convergence of the expansion. 


From the distribution of characteristic values as found in section 
VII, for both the cases of regular and irregular conditions there dis- 
cussed, it becomes readily apparent that given any constant N sufh- 
ciently large it is always possible to choose a circle C whose center 
is at \ = 0, and whose radius both exceeds the value N and is such that 
the distance from any point of C to a characteristic value is greater 
than some fixed constant § > 0. 

We shall consider the convergence of the contour integral in formula 
(95) as the contour of integration is taken successively as a larger and 
larger circle of the type C above. With the size of the circle the num- 
ber of characteristic values which it includes and hence the number of 
terms of the series which are summed by the integration may be in- 
creased indefinitely, the limit of the integral for |\| = ©, being, if it 
exists, the sum of the corresponding series. 

Let us recall the hypotheses already made concerning the functions 


Y:(z). We have 





(i) vi(x) = 0 

(96) 4 (ii) yi(z) continuous a S x S b (page 72) 
(iii) ifn22 yx) + y;,(x) for += 7 3 
(iv) ifn 22 arg{yi(x) — y;(x)} = hi; (page 83) 

To these we shall add 

(97) arg yi(x) = H;(aconstant), 72= 1, 2,...n. 


Concerning the vector to be expanded we shall assume — 


that the elements f;(x), i = 1, 2,...n, consist in the interval 
(98) a <2 3), of only a finite number of pieces each of which is 
real, continuous, and has a continuous derivative. 


It will readily be seen from (96) and (97) that we have restricted 
each +;(x) to vary along a ray in the complex plane. Moreover, if 
n = 2 the dependence of ;(x) upon 2 must be such that the slope of 
the line joining any two of these points remains constant. This 
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means that every n sided polygon with vertices at the points y;(2o), 
7 =1, 2,...n, a S aS 5, at most expands or contracts about A = 0 
as 2 1s allowed to vary. 

If the lines along which 7;(x), 7 = 1, 2,...n, vary are all distinct, 
the conditions of the system are regular provided W, and W, are 
suitably chosen. If, however, one or more sides of any of the polygons 
mentioned lie on a line through the point A = 0, we have the irregular 
case discussed in section VII. In particular all the sides may lie on 
such a line, (for instance the quantities y;(2) may be all real) and it is 
upon this configuration that the irregular case in question bases its 
chief claim for interest. Since the condition (96 iv) is automatically 
fulfilled in this case the functional dependence of y;(2) upon 2 is far 
less restricted than when the points y;(2) form the vertices of actual 
polygons. 

Substituting in formula (86) the value G(x, t) = = Y(x) Z(t) we 
have 

G(a,t,d) = Y@x){ = 31+ 3AtW.Y(@) — W.¥(b)}} 20, 
the upper sign holding for t < x and the lower sign fort > x. It is 
desirable in the following work to express G(x, t, \) in a somewhat 


different form. 
Upon setting 


+31+347{ WY@ — WeY(6)} = 6) + U, 
multiplication by A yields 
W.~Y(a) +4 WiY(b) +3 WY (a) —5 WiY(b) = ; 

{W.Y(a) + Wo¥()} (6%) +4 U, 
which, in view of the relation (5;;) aa (5;;) = J, reduces to 
WY (a) (653) — Wo¥(b) (6;;) = AU. 
Hence 
U = A7W.Y(a) (63) — A Wo Y(b) (6;,). 


Now setting 
—37+4A-{W.Y(a) — WiY(b)} = — (2) + V, 








RBOUNDARY PROBLEMS AND DEVELOPMENTS. 111 


it is found by precisely the same method that V = U. In consequence 
we have for the Green’s function 


(553) 
(99) G(a,t,d) = Y(z) or 7 +47 W.Y(a) G5) 
~(8) 


— A1W.Y(b) (63 Z(t), 


where the upper form is to be chosen when ¢t < 2, and the lower one 
fort > a. 

By means of this formula G(2, t, 4) may be explicitly represented by 
choosing as the solutions Y(x) and Z(t) a pair which are analytic in \ 
and have the forms 


¥(zx) cal (ariorayo 55; " eae | 


Z(t) = (10-80 6:;+ hae | " 


obtained in sections V and VI. It should be observed that G(z, t, r) 
is unique (see page 70), despite the fact that if n = 2 the choice of 
Y(x) and Z(t) thus determined upon changes from any one to any 
other of the sectors within which no quantities R{A{y.(x) — 7;(x)}} 
change sign, and despite the fact that the values of 5; and 53; change 
from any one to any other sector ¢,,. - 

From (95) we have, upon denoting by S»(x)- the sum of those terms 
of the formal development of F(x)- which correspond to the charac- 
teristic values enclosed by the circle C in question, 


b 
(100) Sm(a)- = — ff G(a, t, ¥) R(t) F(t)-dt dr. 


CGC «a 
Substituting in this the value of G(z, t, \) as found above in (99) it is 
4 ' ° 
seen that S,,(x)-= } S“(x)-, where the quantities S(2x)- are given 


t=1 


by the relations 
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1 aie. , , * 
SO (2). = L ve f if Y (x) (0;;)Z(t)R(t)F(t)-dt dd 


Cyuv @ 


b 
Sm(e)-= & = J J Y (x)(65;)Z()R()F(t)-dt dd 


Cur * 


(101) : 
3 ms aa 7 +“ 
S?(z)-= L a f Y(x)A J WY (a)(6;;)Z(t)R(t)F(t)-dtdd 


Cuv 








b 
1 : 
S® (x). = L ae i Y(x)A- | WY(b) (6;,)Z() RF (t)- dtdr, 


\ Cur 


C’,, denoting any arc of circle C which lies within the sector ¢,, and, if 
n = 2, upon which none of the quantities R {\{y;:(x) — y;(x)} change 
sign. If C,, abuts upon a ray bounding ¢,, it shall either include or 
exclude the end point for which the quantity R {AT,,(b)} or the quantity 
R {dT,(6)} vanishes according as the quantity in question is <0 or 
> 0 within ¢,,. The arc C,, may, therefore, include one, both, or 
neither end point, and since the reasoning is precisely the same in each 
case (the sector can be split if both end points are included) we shall 
consider for the sake of concreteness that it includes one, namely that 
one which lies on the ray R {AT,(b)} = 0. The symbol ~ indicates 


that the sum of the integrals over all arcs composing the circle C is to 
be taken. 

We shall proceed to evaluate each of the integrals above in turn, 
and in the course of this evaluation it will be convenient to refer to the 
facts established by the following lemmas. The notation | g(x) |<M 
will be used to indicate that the function in question is bounded. It 
is not to be understood that the M in any one case represents the same 
constant as in any other, but merely that there exists in each case 
some constant for which the relation is true. 


‘Lemma 1: Given any function ¢;(z, \) which is such that 
(i) |gilz,r)|< MforaszxsB, |r| >N, 


(ii) lim | g:(a, \) | = O uniformly, for « < x < B — x, 
A= 0 


(102) | 
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where x is an arbitrarily small positive constant, then 
. 
lim gi(x,r\) dx = 0. 


[| =co 
Qa 


Proof: Under the hypotheses (102) we have 


B—x 
lim gi(a, A) dx = 0 uniformly, 
|A| =o 
8 
and | J ole ) dx |< Mx, 
B-x 
while it follows from these relations and the relation 
8 B—x 8 
| ff ele d) dx | S| ff este d) dx| + | J gilz,d) dz |, 
a a B-x 
that 


8 
| f exlz,» dx | < 2Mx 
for | \| sufficiently large. Since x is arbitrary, however, this means 


that 
B 


lim gi(x, 4) dx = 0. Q. E. D. 


|| = 00 | Z 
Lemma 2: Given any function ¢g2(A) which is such that 


(i) | go(A) | < M for 6 S argd S O3,/A| > N 
(103) “ lim | g2(A) | = 0 uniformly, for@, S arg S 0, - 6,, 
(ii) [A] = 00 
where @, is a constant arbitrarily small but positive, then 
. 1 dr 
ee ea 
Di-® Qaiy oe 
Ce being that arc of the circle | \ | = p which lies in the sector bounded 
by arg \ = 6, and arg A = Og. 
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Proof: Writing \ = pe” we have 


y2(d) = g2(pe”) = 2(6, p). 
Hence 


op 
1 dy ] 
in J wii tim 2 J vx, 0) a, 
Xr p=O te 


(Almoo 272 Rc 


and since the limit of the integral on the right is zero by lemma 1, 


1 
es | ee ex(n) = = 0 Q.E.D. 


IM=co 27 is 


Lemma 3: Given any function ¢3(z, 4) which is such that 
(i) | es(x,d)| < M fora < x S B, 0, S argd S 6g,|d| > N, 


(104) ¢ (ji) tim ¢3(2,d)|= 0 uniformly fora S<2SB—x, 
{= 00 


6. Sargr\ S Os—9,, 





\ 


o 
then I(X) = f ¢3(x, X) dz is a function of the type ge(A) of lemma 2. 


Proof: When X is confined to the sector 0, S arg X S 43 — 4,, 


_ I(A) = 0 uniformly by lemma 1, while for X in the larger 
Al = 00 
sector 6, = argA S 9% 


» 
110) |< f ade = M@ - a). 


Hence I(X) is of type g2(A) by definition. Q. E. D. 

As a matter of convenience we shall use hereafter the symbols 
$1, $2, and ¢3, to designate any function of \ and other arguments 
which satisfies conditions similar to (102), (103) or (104) respectively. 
As heretofore functions which approach the limit zero uniformly as 
|X | increases indefinitely will be denoted by €. Let us return now to 
the direct evaluation of S,, (2)-. 
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I. S@(a)-. 


Writing 





I= — f ¥(@) Gi) ZO RO FO-a 
we have from (98) 
1 
S(z)- = — f svar. 
271 Cc 
From this it is immediately seen that 
(105) SY (a)-= 0. 


Consequently it is necessary to consider further only the case z+ a. 
We have 


Ji = ( : 3 OTOL 15 41 /Npie(x)]}opme Or om 


k,m,p,q=1 
LSmp + 1/AYmp(t)]}dpgve(t) fa(t) it) : 


le. 


(106) Ji-= — ( f MPO RMA Sirvi(t) fill) it) oo 


a 


V/n( XL eM TEMPO} FE RO FF volt) fp(t) [bap vie(z) + 
k,p=1 


Dik Vep(t)] it . 


For \ on any arc C,,, however, each term under the sign of summation 
in the last matrix on the right of this expression is seen to be of the 
type v(t, \). Their sum therefore is of the same type, and by lemma 
2 the integral is of type g2(A). Consequently we have, upon integrat- 
ing by parts the elements in the first matrix on the right of (106),*° 


sinbies ‘(- bis — fi — 0) +e iF (a +0) + 


° (x)—-T. d r)-B.( (ga) 
{Pwr © \ ¢ yy B(-BA) 


15 Recall that r;(a) = B;(a) = 0. 
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Accordingly it is seen that J,- can be considered as a sum of ma- 
trices of which the first three are discussed as follows. 
For the first matrix we have directly 


l * ] ea ; 
x (55; f:(x —0)) = x (5;;) F(z — Q)-. 
The second matrix 


] * (zx Az 
ca : (5, eril )+ BC ) fi(a)) 


is one each of whose elements is seen either to vanish or to be of type 
g2(A)/A on are C,,. Hence the entire matrix may be represented 


by (y2)/ d. 
The third matrix 


] ( d 


is likewise one each of whose elements either vanishes or is, on arc 
C,,, of type g2/ by lemma 3. This matrix is, therefore, also of the 
type (¢2)/A, and inasmuch as a sum of matrices of type (¢2) is again 
a@ matrix (g2) we have 


pee : {(3%.) F(z — 0): + (¢s)} 


The integration of J,- over the arc C,, yields, therefore, by lemma 2 
1 rf * dy Ww rina 

— (uv) a ; —=-— 5;* F(a — 0)-+ (6), 

2771 J wig Mbdtincdtts: Windle rd 2a sialiidie ) 


where w,, is the angle subtended by C,, at \ = 0. 
A similar expression results from integration over each arc similar 
in character to C,,. Hence we have 


SP (a)- = —” (6) F(x — 0)-+ (6). 
Cc 2r 
Let us consider in detail the sum 


(107) xr == ef”). 
C 2r 
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Corresponding to each arc C,,, there exists another arc ';; which is the 
reflection of C,, in the point \ = 0. Since C7; also subtends the angle 
w,, at the point \ = 0 we can write the sum (107) equally well in the 
form 


(108) Xr =” (5 + 5G), 
10 Q2r J 


the summation covering now the arc of only one (any) half of circle C. 
* x 
But it will readily be seen that 4” = 6%”, for quantities which have 


a positive real part in o,, have naturally a negative real part in o;;. 
* bei 

Hence eat +- ik = 6,;, and the sum (108) i.e. (107) reduces to 

XI =31. 

$C 27 
It follows, therefore, that 

Wuy (u») 
(109) YS OF") = 71. 
Cc 2T 

In consequence we have 


Sin (x): = 3 F(x — 0)-+ ©), 
i.e. 
(110) lim S(x)-= 4 F(x — 0)- 


m= 00 


II. S®(z)-. 





The treatment of this expression is parallel to that of S\)(x)- and is 
as follows. 


Writing 
b 
Ja = J Vw) OZ © RO FO-< 
] 
we have from (101) §@(x)-= = 
2m % 
Hence 


(111) §@(b)- = O, 
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and it is necessary to consider further only the case x + 6. We have 


Je ae (f > oly \2)+ By (2) 5; ul. Soc wate) os oor m-B,, = bn 4. 


k,m,p,q=1 


Ym(t)] 
r 


b 
Pi (J MEMO HOO FF yilt) fill) it + 


x 





pcre) fa(t) i) 


i = op ve 
" ( & onl Pr, (2) r,(¢) }+B,( ) —_ Y(t) fp(t) [Sip vin (2) 
2 p= 


+ bik Vip(t)] it) 


Each integrand in the last matrix of this equation is of the type 
gi(t,) and accordingly the matrix is, by lemma 1 of type (€)/d. 
Integrating by parts the elements of the first matrix on the right we 
have 


l ** 
alia (3 — MMOMO +B (2B £5 — 0) + file +0) + 
b 
_rig; @ (e) 
Je 7 fe + “' 


Therefore, J2- can also be expressed as a sum of matrices, the first 
and third of which, namely 


1 (— af ater nin omc-90.) 


and 


b 
L (oe (° xgraer-rioy 4 4 > z)- : 
“ (3: | ert P(z)-T;( } = { f(t) Bi ) B,C) at), 


are of the type (e)/A on arc C,,, for each element either vanishes or 
approaches the limit zero as |A| = ©. The second matrix of the 
sum is directly 
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~ (6% fle +0) = = (i) F@ +0)-. 
Hence 
Be ; {(8%*) F(z +0):+ 0}, 


which integrated over C,, gives 

1 oe 
sdf (5%) F(a +0)-+ @\ > aH (5) F(x + 0)-+ (6). 
2m Cu» 

By familiar reasoning this leads to the equation 


SH (x): = 2 =~ (3) F(x +0)-+ (6, 


‘T 


and inasmuch as 


(112) is by (a) = 1], 





bo 
+ 


as may be shown by applying again the argument by means of which 
relation (109) was established, we have 
Sin (x): = 3 Fe + 0)-+ (©), 
1.€. 
(113) lim S®(x)-= 4 F(x +0)- 


m= Oo 


WI. S®(z)-. 





Following the procedure used in the discussion of the preceding 
expressions, let J3- be defined by the equation 


J3-= — Y(z) A i WY (a) (6;;) Z(t) R(t) F(t) -dt. 
l 
Then SY (x): = ay Pe 
2rt % 


se 16 
(C8) 





16 See note on page 104. 
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and substituting for D(A) its value as given by equation (56) we have 
; n eK 
ar el *Te(b)+ By (b)} 5”? Dji(d) 


k=l 
A- Pe kt+?t 











n *K 

ar ol Te(0)+By)} BH” Dy) 

| a 

’ 
os oz 

ae o7 {AT {O)+B,6) } 5(uv) D;d) 


D(a) 








\ 


where D(a) is given for the case of regular conditions by formula (57) 
and for the case of irregular conditions discussed in section VII by 
either formula (66) or formula (68). A glance at these formulas 
shows, however, that in each case the expression for D(a) reduces, 
under the conditions for which the expression is valid to the form 


D(s) = W™ + g2, on are C,,. 


Now every point of a circle C, and hence in particular of an are C,, 
as at a distance which exceeds a fixed quantity from any characteristic 
value. However, by analogy with formula (59), for any point of C,, 
under regular conditions for which the function D is sufficiently 
small we have 

“es — yw) 
A= r,(b) } - B,(b) + log OT + e+ 2prif. 

In other words the point in question will necessarily lie near one 
of the characteristic values. 

This stands in contradiction with the fundamental property of the 
circle C. 

Inasmuch as a similar relation holds on every arc of type C,,, 
it follows that for every point of all circles C, D(A) exceeds, under 
regular conditions, some fixed positive constant 41, i.e. 


(114) | D(A) |2 6: > 0 for d on any circle C. 


While we have thus proved the relation (114) only for the case of regu- 
lar conditions it can be shown to hold equally well under the type of 
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irregular conditions considered in section V II.” In either case, there- 
fore, A“ is uniformly bounded on all circles of type C. 


Proceeding, we have 
f . 
Y(2)A7= > orl y(@)+ By) bat me eae t _{AT,0)+B,0)3 8 Die) 


ft DO) 


f ‘ 


eM T(2)-Ty0) 85 J+ Bi@-B,O) 5. tol 














D(a) 
j 
eu > en af Py(2)-1, (015.4 }+By(2)-B,(0)34* [esn(x)] Djx(d) 
j 





Since the last matrix of this expression is of type (€) we have 


Y(x)A1= (a 1 (2)-14(0)8 54 }+B,(2)-B,(0)3;, } 2a) + €E ), 
D(a) 


or, more briefly, 
(115) Y(x)A-1= (el I',(z)-1,(b)§ 45 1+ B,(2)-B,(b)8;; ris(d)), 
where — [| re) | < M for |A| > N. 
Inasmuch as D(A) was seen to be of the form 
Dr) = W™ + o 
for \ on arc C,,, we have for any X on this arc 


Wy) be W (uv) 
(743) = (; a + . = (“ 7 2) + (¢2), 





or, denoting by 2” the matrix whose determinant is VW”, 


~1 
(116) (ri;) = 2” + (¢»). 





17 Cf. Wilder, loc. cit. p. 422. 
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We en further, 


J WY (a) (83;) Z(t) R(t) F(t)-dt = 

















Sy i s-lemtel} ap Mf g(0-By | oe 
h,k,m,p,q=1 ; ita 
ott 
Yeaill - V(t) f(t) dt) - 
b 
E wlan f PVM 409 fy at | + 
h=-1 
; ( : : . 
" 2 wih bak , e TKO FEO y(t) fo(t) [onWep (t) + Sep gre] at 
»k, p= 
' 
f 
1 
= : > ws Os — eNO) By) ¢ (b—0O) +f,(a +0) + 
h=1 





(6) 
2 





b 
fe or, 4 ) fa(t) HO a 


j 


the last form representing as heretofore the result of an integration 


by parts. The multiplication of this matrix by the matrix (115) 
yields the result 


d 
nm 


— | ~ * es 
re es s = el r';(z) P;(d)6.. , +8; (2) Bia se7(d) ws wie Derk 
=] 





\ 


b 
fa be 0) mA EMO B®) 6 (5, — 0) re f emi S {f,(t)e 2h } dt 
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the principal matrix on the right again breaking up into three compo- 
nent matrices. Of these the first, 


n x ko 
1 ( 3 Mf P(x)-T;(6)6 ii }+B(z)— BAS 5 rix(d) we orn ert a(b)— B00) fulb—0)), 
 \k,h=1 
and the third, 


4 

on " +e + 

poets... > er Pj(x)-T (0) 555 }+B,(2)-B()5;; 
N | & het 





b ) 


Tin(d) WP Onn fe anya Lene (t) }dt 


‘“ 





j 


are for all values of a clearly of the type (€)/A. Only the second 
matrix of the sum, namely 


a aad & (2)-9.(b)a. 14-B.(2)- B.(b) 820 
J3(x)- = — ( | a cael i a acetate ras(d)wiQ bis fa(0-+0)) 
kh=l 
remains, therefore, to be considered. It depends for its character 
upon the value of x. 
For « + a, x + b, the matrix is clearly of type (g2)/A. For x = a, 
on the other hand, the exponential factor common to the elements of 


KK 
th -f{Xr- * 
the i" row reduces to e7(*7i+F: }sit — 55. + € so that 


J3(a)-= —( 5 > 56 Tsx(A) WY) Sin fx(a +0) + c). 


k,h=1 
The fact that 
J3(b)- = — (5 Y 555 76e(A) wR Sha faa +0) + e:) , 
k,h=1 


may be established in similar manner by use of the relation 
* 
el TiO+B(0)}55; — 5" 4. oy. 


Substituting for (7;;) its value as given by formula (115) it is readily 
seen, therefore, that on are C,, 
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= «ui * —1 ** 

Js(a)- = — (6) oO” W, (6) F(a + 0)-+ (¢2)} 
i will 1 -- * 

Js(b)-= — {(o%”) Oo W.(5Y”) F(a +0)-+ (2)}. 


It follows from this that 





ee 1 fe when x + a, x + b. 
C 271 Dus Xr 


1 rs -l *K AN 
=> ae 2 {—(84”) & W. (6%) F(a +0)-+ (g2)}— 
(3) c 2 C,, d 
Sin (x)- ' 
when x = a 


7 RE. ex 
= 51 J 1-0) 2 1.0) ra+0)-+ Clee 





when x = J, 


or, upon defining the matrices K; and K;; by the equations 


ae tir rile clan\ wa tale 
Ri= 5 }— OPO w. OP) 


(117) . ; 
= W uy pod ay ad 
| K; = x { — 2 of) 0 wr. aff, 
= 0 when zta,x+ b. 
(118) lim S(x)- ( = K3F(a+0)- whenz = a. 


= K;F(a+0)- whenz = b. 


It should be observed that the values of K3 and K; depend only 
upon the differential system whose characteristic functions form the 
terms of the expansion in question and are in particular entirely inde- 
pendent of the vector F(z): 


IV. S*(2)- 





Proceeding precisely as in the case of Sin (x)- we have upon defining 
J4- by the relation 
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; b 
Je= (2) 44+ f Wi@) Cy) ZO RO FO-de, 


1 
SA (x). = 2 f sear 


Further 


J wre Gi) 20 RO FO-a = 


b . 
nr 
Ew aia f AVVO A (0 a 
=|] 


a 
\ 4 


ar b 
Eu she, | MOO 1,4) £500 adap O 
a 


h,k, p=1 








~~ i 





. ’ 


+ Sip onx(b)] at 





L wy in| — fx(b — 0) + PTAO*FOF, (a +0) + 
=] 


a i= 





\ 


| 
ferment { ¢By(b)- By f;,(t) in| he ee 





a 
j 


and multiplying this matrix by the matrix Y(x) A“, as given by 
formula (115) it follows that 


a 


“ “ 
| EAN NOBE HABE y(n) wba | —fa(0—0) + 


\ 


Je= 





, 


b 
d & (v2) 
“Aeaitsio fil a +0) he J er{T,© I, (t) } eo f,(t) ePr') By) dt sa 





4 
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The second matrix of the sum on the right, namely 


1 * +* = ek 
: ( Y AMO MMs }+B(2)-BO)5i57 (rn) ws ner et af, (a+0)), 
k, h=l 


is seen to be of the type (g2)/A, whereas the third, i.e. 
| ~ bt tk 
3s el P;(z)-T;(6)5;; }+ B,(z)-B(b)5;; 


1 
X | kha! 


\ . ; 


d 
raul atx J oH) =, { fa(t) e-PH} at 


a 





+ 





4 


is also of this type, as is apparent from the fact that by lemma 3, each 
of the integrals with respect to t which actually occurs in any element 
of the matrix is of type g2, while the remaining factors are bounded. 

The remaining matrix of importance in the expression for J¢-, 





namely 
_ —|] | “ Od *% * 
Ja(x)- = — s el P-TiO3 4; }+ B@)- B55 7 (nw 55, f,(b—0) | » 
k,h=1 
{ 4 





again depends for its character upon the value of x, being of the type 
(g2)/X when and only when xa and «+b. By reasoning now 
familiar the results 


_ sar 4 1 8 7 * 
Js(a)-= - "i ( x 65; Té(A) WEP Onn dn(d — 0) + e) 
»A~l 


and 


_ Jab): =~ ( L 855 Tin(d) WER Sin tr (b — 0) +.) 


NX \g het 


—1 * l 
— (5;;) QO) Wy (6;;) F(6 —0)-+ (yo), 
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may be deduced, whereupon it follows that 


| 1 dn 
i oy ey. when z + a, x + 8, 
2m ¢,, 


1 f ES nee dd 
sete, | = Fgh (om 1 ab) 10-014 oot 2 
when x = 4a. 

l * 
a ELS ay W »(5%”) F(b ae: 0)-+ (oo) > 


when x = BD. 





Defining the matrices K, and Ky by the equations 


an a3 7 -1 * 
Ki= 2 —— (6) o&™ W, OY”) 
(119) al 
BE §—*" afnod wy 
Ki= 2 (oy) 0” W, apt, 
it follows, therefore, that 


= () when z@ta,x+b) 
(120) lim S“@(2)- { = K,F(6—0)- whenz=a 
tii = K,F(b—0)- whenz= b. 


A summary of the various results as contained in formulas (110), 
(113), (118) and (120) is seen, in virtue of (105) and (111), to yield 


lim S,,(x)-= 4 F(a — 0)-+ 3 F(a +0): when zta,x+ 8, 
lim S,(a)-= 3 F(a +0)-+ K3F(a + 0)-+ Ks F(b—0)-, 


1 F(b — 0)-+ Kz F(a +0)-+ Ka F(6—0)-, 


lim S»(b)- 
In consequence we have the following , 


THEOREM: Given that L is any vector differential system of the 
type 
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Y'(x)-= {A(x)A + B(x)} ¥(z)-, 
W. Y(a) me 5 W, Y(b)- = 0, 


which can be reduced by a change of the dependent variable to a 
system of type (46) for which (a), R(x) and B(x) are continuous to- 
gether with their first derivatives (6), thefunctions 7;(7), 7 = 1,2, ...n, 
satisfy the relations (96) and (97), and (c), the condition (ii) on page 
89 is fulfilled. Then the development in characteristic functions of the 
system L which is associated with any vector F(x): whose elements 
satisfy condition (98), converges to 


4 F(a —0)- +4 F(a+0)- whenzrzta,r+b, 
H,F(a+0)-+ J,.F(6 —0)- when z= a, 
H,F(a+0)-+J,F(6 —0)- whenz = 8, 


the four matrices of constants Ha, Ja, H», and J, being explicitly 
determined by the matrices R(x), A(x), W, and W, as stated. 

















